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Executive Summary 
and 

Conclusions 

PART A; SNOW STUDIES 

Snow events in the American River basin and subsequent runoff condi- 
tions are described for two years: a "typical" yeaf-> 1975; and a record 

drought year, 1977. For the 1975 data, a graph of snowcover plotted 
as a function of Julian day gave a correlation of 0.99+ for the seven 
measurements, while 11 NOAA VHRR measurements gave a correlation of 0.97. 

Thus, the use of Landsat or VHRR data to forecast the date of snowfree 
conditions in the river basin is a distinct possibility. 

The use Of satellite images in conjunction with ground based 
measurements can increase the effectiveness of getting the drought 
"message" across to the public. 

The differing response of various sensors to snowcover measurements, 
previously demonstrated by several workers, should be taken into account 
by snow mappers. For example, 10 VISSR measxirements of snowcover in the American 
River basin are, on the average, 3.5% lower than corresponding values 
from VHRR. 

A comparison of conventional and computer-generated analyses of 
snowcover disclosed that the snow-cover area (SCAj ratio to computer-generated to 
conventional analysis ranged from 0.60 to 0.08. The ratio consistently dropped 
as a function of Julian day and was always least just before complete 
snowmelt. These findings were not unexpected^ as pixel counts invariably 
are more precise than techniques that merely outline ax'eas of snow. 

Further, the smaller the area of snowcover, the greater the error potential; 
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i.e., as the area o£ snowcover approaches zero, the percent error increases 
if the absolute area of measurement in error remains the same. 

An analysis of snowcover extent by the GE Image 100 on a pixel by 
pixel classification, using MSS bands S and 7, resulted in much lower 
snowcover area figures. It was not possible using this technique to 
distinguish metamorphosed or "melting” snow areas, unambiguously. The 

current MSS bands are not optimum for measuring areas of metamorphosed 
or melting snowcover. 

The low values of SCA as measured by coinputer-generated thematic 
mapping indicate the difficulty in properly defining the various themes 
representing SCA as viewed by Landsat. Previous studies have shown that 
conventionally derived NOAA/NESS SCA estimates agree well with ground 
observations. Nevertheless, even if the computer-generated thematic 
system produced snow maps equivalent to conventionally prepared maps, the ; 
computer-generated technique would be more time consuming and costlier at 
this stage. 

Rather than classifying melting and non-melting snow, the thematic 
approach had instead apparently divide.d the melting or metamorphosed 
snow into various degrees of metamorphism: actively melting, early 

melting^ and incipient melting. If this be true, the need for more detailed 
studies of the radiative properties of melting snow becomes more and 
more clear. 

In a study of snow in forests, Landsat MSS when compared with NOAA 

VHRR visible tnagery had a better signal -te-holee ratio and better sensitivity 
Careful digital analysis in the Adirondack Mountains of New York indicate 

that forest type has a significant impact on reflectance. Deciduous and 
conifer areas show increases in radiance as a function of increasing wave- 
length. The converse is true for open areas. As forest areas are critical 
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for snowcover studies, careful controlled ground measurements of radiation 
are needed. 

Snowcovered area was shown to be an important variable in runoff 
prediction by a stepwise multivariate analysis for several American 
River subbasins. Satellite-derived | SCA and day of year provided estimates 
of daily average runoff with a correlation of about 0.90, during seasons 
in. which the recession curve of the hydro graph is fairly smooth. 



PART B: SOIL MOISTURE STUDIES 


There is a real need for areally averaged soil-moisture values 
that reflect the condition of the watershed and its response to rainfall 
more accurately. We attempted to relate field studies of soil moisture 
near Phoenix, Arizona, and Luverne, Minnesota, to Landsat spectral 
reflectance measurements not only with ground truth but also with 
conciirrent gamma -ray aircraft data and cooperatively with NASA PMIS 
(microwave) overflights. 

At Phoenix, soil moisture values computed from the gamma-ray gross- 
count rate of the 18 March 197S data are slightly higher than the ground- 
based measurements. The mean of the ground-sampled moisture values was 
nearly the same as that computed from the photopeak gamma-ray technique. 

At Luveme, the gamma-ray adjusted gross-count -rate method gave 
an average soil moisture value of 28%, with mile-by-raile values ranging 
from 21-36%. Ground-based soil moisture measurements gave an average 
of 28%, with individual values ranging from 23% -36%. Thus the gamma-ray 
adjusted-gross-count technique provided a reliable source of corrobora- 
tion for the areal soil-moisture values. 

Gamma -ray aircraft surveys provided only one value for the entii'e 
Phoenix test site as a result of averaging many data points Whose 
"footprint" was 150 m. The technique is limited by atmospheric contami- 
nants, the need to fly at low (150 m) altitudes , and the need for prior 

■ ... . I ■■ ■ . ' ■ ■ ...... .... 

ground calibration of a given test site. 

Multispectral data can be used to establish a spectral signature 
for specific soil types. It is possible to detect significant changes 
in reflectance and relate these changes to changes in roughness 
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(i.e., tillage), soil moisture, or vegetative status. Identification 
of the causal factor is equivocal. 

A comparison of MSS band 7 raw data (full resolution) over the 
Phoenix site superimposed on a false color near-lR aerial photo taken 
27 hours later revealed that for alfalfa field 97A 0-2 cm soil-moisture values 
were lowest (13.70%), and band. 7 reflectances were lowest (64). For alfalfa field 
97B 0-2 cm soil-moisture values were highest (27.65%) and the band 7 
reflectances were highest (84) . This direct relationsh.ip of soil moisture 
to hand 7 . values held for all alfalfa fields sampled in the study. 

The amount of useful information from visible and near-IR data is 
miniscule when compared to that provided by passive microwave radiometers 
and, potentially, from the therraal-IR. Further efforts to estimate soil 
moisture from remote platforms should place the emphasis on the microwave 
and the thermal IR portion of the spectrum where the most encouraging 
results have been achieved. 

The use of the Geostationary Operational Environmental Satellite 
(GOES) thermal IR system to monitor soil moisture will not be feasible 
mtil the spatial resolution is improved from the current 8 km, and 
unless computer compatible tapes are made available. 


INTRODUCTION 


This report was prepared as a result of a grant to the National 
Oceanic and Atmospheric Administration/National Environmental Satellite 
Service (Investigation No. 23170, Contract No. NAS-53991A) . It is pre- 
pared in two sections: Part A, which deals with snow studies, and Part 

B, which deals with soil moisture studies. The purpose of this report 
is to evaluate our ability to extract useful information from Landsat-2 
data, and to apply it to current needs of the National Oceanic and-. 
Atmospheric Administration's hydrologists. As can be seen throughout 
this report, the authors have continually tried to examine Landsat data 
and compare it »vith the current operational NOAA series of polar-orbiting 
satellites and the geostationary series (GOES). 

The authorship of the report is a team effort, but nevertheless, 
certain sections are primarily authored by individuals. In Part A, 
section one and two were written by Wiesnet, sections three, four, and five 
were written by McGinnis. Sections one and two of Part B were written I 
by McGinnis and section three by Matson, Each has i 

reviewed the others' contributions and the ranking of each author represents 
the amount of time (work) spent on the entire project. Please note that 
the references are found at the end of each section for the readers' 
convenience. 

Lan,dsat-2 was in no way different from Landsat-1 with respect to 
sensors. We have stressed the Multispectral Scanner or MSS data. Our 
results differ from the Landsat-l report (Wiesnet, et ai. ; 1975) in the 
fact that we have narrower objectives, i.e., snow and soil moisture, and 
in the increased use of computer compatible tape data (CCTfe), rather than 


the use of photointerpretive techniques. 

We lament our inability to obtain aircraft/satellite coordination on 
our experiments, but all concerned must realize that the acquisition of 
NASA aircraft data precisely in step with the satellite- is indeed very 
rare. The limitations of project priorities, weather, and equipment 
failure make coordination of aircraft overflights a nearly impossible 
task. We were able to get aircraft data from gamma-ray overflignts 
;performed by EG§G, Inc. under contract, and these data have been incluaed 
in Part B. I 

The access to NOAA imagery and image processing systems has been 
of great help to us, and we have freely used and relied on NOAA data, 
computer facilities, and equipment. Whether this report answers more 
questions than it raises is for the reader to determine. It certainly 
offers some interesting conclusions, but it is also clear that many 
unanswered questions remain. 
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INTRODUCTION 


Snow, being the supplier of water for human, agricultural, and 
industrial consumption, is a vital water resource in arid areas of the 
world, A complex reservoir system, such as that in California, effec- 
tively transports snow runoff to meet the demands of cities and farms. 
However efficient these systems, better snow/runoff forecasts are 
always needed. The snow studies presented in this report attempt to 
evaluate the role satellites, principally Landsat, can provide in 
improving these snow/runoff forecasts. 


Four chapters discuss various snow parsimeters and problems encountered 
with snowpacks in mountains. Section 2 evaluates several years of snow 
extent mapping in the American River basin — a review and update of ouir 
Landsat-1 study (Wiesnet, et al., 1975) . Snow- extent mapping from satel- 
lite imagery has successfully evolved into an operational mode with 30 
basins in the itfestem United States currently evaluated by NESS through- 
out the melt season. Snowmelt in the American River basin is reevaluated 
in Section 3 and updates work begun in our first Landsat study. The major 
emphasis in this study incorporates analyses performed on the G.E. Image 
100 Interactive Multispectral Image Analysis System. A summary of the 
work done by Eschner, et al. (1977) regarding snow in forest areas 
performed as an adjunct to this study constitutes Section 4. Multiple 
regression analyses for three subbasins in the American River basin 
relating snowmelt runoff to various meteorological and satellite 
measurements are investigated in Section 5. 

This portion of the report represents another step toward application 
of Earth-orbiting satellite data to the real day-to-day problems of water 
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resource studies, research, and management. As a research paper it does 
not attempt to be comprehensive, rather it attempts to demonstrate by 
simple exauaples, that new ways of assessing or monitoring the snowpack will 
come only after we fully understand the type, quantity, and quality of 
information inferable from satellite orbit. 
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Snow Extent Mapping 


The capability of Landsat MSS and NO(\A VHRR to measure the areal 
extent of snow in river basins on a repetitive basis has been amply 
demonstrated. (Barnes and Bowley, 1968; Wiesnet, 1974; Wiesnet et al., 
1974; Rango, 1975; Meier, 1975; Schneider, et al., 1976). This section 
will report in some detail on the snow extent in the American River basin 
which lies on the western slopes of the Sierra Nevada in California just 
west of Lake Tahoe (Fig. 1). In our original proposal, we had 
planned to include the Feather River basin as well, but this latter area 
was also included in a snow study by A.J. Brown, California Department 
of Water Resources, under the Advanced Systems Verification Test 
(ASVT) program of NASA (Rango, 1975, p. 0 ), Therefore, it was deleted 
from our plans as duplicative. 

The fertile Sacramento and San Joaquin valleys of California depend 
heavily on the Sierra Nevada snowpack to provide water for agricultural 
use as well as for municipal and industrial use. Nearly 50% of the 
average annual runoff takes place during the April-Juiy snowmelt. Tlie 
California State Cooperative Snow Survey Program began in 1929 to 
measure snow and to forecast hydrologic runoff. Water management in 
California is as sophisticated as anyxvhere in the world. The State has 
been quick to recognize the potential of satellite data collection and 
satellite snow observations (Brown, 1975) . 

A. The American River basin 

1. Description . This rugged mountain. basin lies between the 

O 

city of Sacramento and Lake Tahoe. It covers 5,601 km . The Central 
Sierra Snow Laboratory lies just outside the basin boundary (Figure 2) . 
Much of the basin receives more than 152 cm (60 inches) of precipitation 
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largely in the form of snow. In the spring and summer, clear skies 
predominate and facilitate satellite observations. Runoff is about 102 cm 
(40 inches) per year, and reaches a peak during the April-June period. 

2. The 1975 water year . In the American River basin the snow- 
pack built up slowly in November through January, but heavy snows fell 
in early February, changing the outlook for water supply from near drought 
conditions to a higher- than-average year. Figure 3 shows the seasonal 
distribution of snowfall. The snow=?extenu build-up recorded by NOAA 
A/HRR images (figure 4) shows a dramatic rise during the first week of 
February from 22% to 62% in the snow-covered area (SCA.) . The rapid 
decline thereafter was due chiefly to lovrer level melting. Once into 
Spril, the melting of the snowpack assumes a fairly stable rate. Snow- 
fall and snowstorms decrease in amount and number. 

The Landsat MSS record of melt (figure 5) appears much smoother 
than the NOAA VHRR, because of fewer observations, but is at about the 
same rate. 

Comparison with the 3-day running mean of the NOAA VHRR-derived SCA 
data shows even closer similarities (figure 6) „ Comparison with the 
GOES values (figure 7) shows that GOES consistently detects less snow 
than Landsat; this phenomenon is attributed to the different spectral 
response of the two sensors; GOES peaking farther in the near-IR receives 
less visible radiation (Schneider and McGinnis, 1977). 

3* The 1976 water year . This was the ghird driest year of the 
century in California. Only about 76 cm (30 inches) of precipitation 
fell over the American River basin as an offshore high-pressure system 
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blocked the movement of cyclonic storms into California. This figure 
is only 50% of the mean. Although the Blue Canyon weather station 
received 80% of normal precipitation in February, by 1 April the water 
equivalent of the snowpack in the American River was down to 32% of the 
mean. For the entire water year, the unimpaired runoff of the American 
River amounted to a scant 31% of the mean runoff. Table 1 shows the 
comparison of the water supply forecast with the observed flow. The 
difficulty of forecasting during extreme events such as prolonged drought 
is clearly evident here from the figures in the table. 

Figure 8 shows the water equivalent of the snow at a high- j 

altitude station. Note the steep dropoff to zero in mid-May, and the 
average dropoff, which reaches zero on 20 June. 

4. The 1977 water year . The second drought year in a row 
plagued California. Runoff from the North Fork of the American River at 
North Fork Dam — ^an index station--was only 8% of median. The December 
monthly mean dischargeof 33.1 cfs and the daily mean of 29 cfs on 
28 December were the lowest in 67 years of record. The drought was 
reportedly the worst in almost 100 years. In January the North Fork Dam 
station's monthly mean rose to 63.2 cfs which was still the lowest of 
record for the month in 36 years. Folsom Reservoir at month end was 
only 30% of normal maximum. The snowpack was far belovi/ the mean. 

In February the snowpack increased slightly but strearaf low remained 
below normal for the fourth consecutive month. The North Fork Dam gaging 
station recorded a monthly mean discharge of 77.2 cfs ( 6 ^ of median) , the 
lowest for the month in 36 years of record (see figure 9). Folsom j 
Reservoir dropped to 27% of normal maximum. 

In March the North Fork Dam gage recorded new low- flow records again. 
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Table_J_» Observed and Forecast For the American River inflow 
to the Folsom Reservoir (in acre-feet x 1000) (based on data from 
California Department of Water Resources) 

Forecast 


Date 

Period 

Flow (efs) 

Observed Flow 

Forecast 

Feb. 1 

Apr-Jul 

420 

312 

+35% 


Water yr. 

1,120 

785 

+43% 

Mar 1 

Apr-Jul 

485 

312 

+55% 


Water yr. 

990 

785 

+26% 

Apr 1 

Apr-Jul 

300 

312 

- 4% 


Water yr. 

755 

785 

-10% 

May 1 

Apr-Jul 

245 

312 

-21% 


Water yr. 

700 

785 

-11% 
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WATER CONTENT, inches 









The monthly mean of 109 cfs (-9% of median) and the daily mean of 73 cfs 
on the 6th were the lowest for the month in 36 years of record. Folsom 
rose slightly to 29% of normal maximum. 

In April new monthly lows were recorded. The month- end snowpack was 
far below normal (figure 8). Folsom reservoir was 30% of normal maximum. 
In May, the sixth consecutive monthly low flow record was set at North 
Fork Dam. Folsom Reservoir held at 30%. June was the ninth consecutive 
month of below-normal streamflow at North Fork Dam. Folsom dropped to 26% 
of normal maximum. As of 31 July, cumulative runoff for the first 10 
months of the water year was only 12% of median! Folsom Reservoir con- 
tinued to drop to 24% of normal maximum. 

In August, the flow at North Fork Dam was 6 cfs less than the 
previous record low August mean (1931). The cumulative runoff for the 
first months of the water year was 12% of median. Folsom Reservoir was 
only 17% of normal maximum. 

The monthly discharge of 14.7 cfs was the lowest of record for 
September, and the annual mean discharge was 88.4 cfs for the 1977 water 
year“-the lowest of record. 

5. Landsat MSS coverage and NOAA VHRR coverage . In 1975, 
during the snowmelt season, seven MSS measurements of basin snow extent 
were made (see Table 2). Only one pass was obscured by clouds. During 
the same period, 11 measurements were made using VHRR images, and 10 
measurements were taken drom VISSR (GOES) images. 

Figure 5 shows the snow-covered area depletion curve prepared 
from the Landsat MSS band 5 images. A simple regression of snowcover 
plotted as a function of Julian day yielded a parabolic correlation 
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TABLE 2 


Lands at -derived snow 
1975. 

extent measurements 

for the American River Basin, 

Satellite 

Date 

Percent Snowcover 

Lands at -1 

4/29/75 

49 

Lands at -2 

5/8/75 

45 

Lands at- 1 

5/17/75 

40 

Landsat-1 

6/4/75 

24 

Landsat-2 

6/13/75 

16 

Landsat-1 

6/22/75 

12 

Landsat-2 

7/1/75 

7 
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coefficient (r) of 0.99+ for the seven MSS measurements. Using 11 
points measured from the VHRR, the same plot yielded an r of 0.97. 

Thus the use of Landsat snowcover data to forecast the date of snow-free 
conditions in the river basin is a distinct possibility. The major 
deterrent is the slow delivery of Landsat data to the user. The more 
timely delivery of NOAA VHRR-based SCA maps makes them a preferred source, at 

present. 

6. GOES VISSR coverage . The NOM satellite. GOES (Geostationary Operational 
Environmental Satellite), is useful for mapping snow by means of its 
Visible Infrared Spin Scan Radiometer (VISSR). The satellite system 
consists of two satellites over the equator, one at 75°W and the other 
at 135°W. The VISSR has a visible band (0.55-0. 75iim) and a thermal IR band 
(10. 5-12. 5pra). Ground resolution of the visible band is roughly 
equivalent to the VHRR (1 km) . The resolution decreases as the "look" 
angle increases; thus the far north (above 55°N) is not a suitable target 
for snow mapping. The oblique viewing angle also creates some panoramic 
distortion in the GOES images. Nevertheless the VISSR is capable of 
producing an image of the same area every 30 minutes, a feat no other 
satellite sensor performs. The advantages of VISSR data then are that 
it greatly increases the probability of securing a cloud-free image 
of a given river basin, and it permits multiple looks on a specific day, 
thereby revealing the effect of sun angle on snowcover detection and 
reflection, 

A plot of 10 VISSR-derived measurements of snowcover are shown on 
figure 7. Snow-surface reflectance using the broadband VISSR is 

consistently less bright and less extensive owing to the near-IR shift 

of the VISSR spectral response curve (Schneider and McGinnis, 1977). The _ ‘ 
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reflectance of snow is greatly diminished in the near-IR as demonstrated 
by Strong et al. (1971), Wiesnet et al. (1974), and O'Brien and 
Munis (1975) . Figure 10 demonstrates that the GOES VISSR consistently 
records less snow than Landsat MSS Band 5 and NOAA VHRR visible band. 

In 10 measurements of the American River basin the VISSR read an 
average of 3.5% lower than corresponding VHRR values for the same day. 
The differences ranged from -1.5% to -6.5%. 
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Snotvmelt in the American River Basin, California 
Introduction 

Several studies (Wiesnet, et al., 1975, Thomas, et al., 1978) 
hiave attempted to derive various snowpack parameters using the multispectral 
approach available ivith Landsat MSS data. Early techniques using raw, 
uncorrected MSS reflectance counts (Wiesnet, et al., 1975) have been 
superseded by more sophisticated analyses that use computer programs 
(Thomas, et al., 1978). Results thus far are mostly preliminary, although 
measurement of snow extent has become routine and operational. Informa- 
tion concerning the density, moisture content, and the surface condition 
of the snowpack derived from MSS data is still speculative. In our attempt 
to extract information on snowpack variables from Landsat MSS computer 
compatible tapes (CCT's) of snow-covered scenes, the GE (Generai Electric)| 
Image 100 Interactive Multispectral Image Analy,sis .System was_ utilized, j 
Description of the G.E. Image 100 System 

The G.E. Image 100 Interactive Multispectral Image Analysis System 
(henceforth called the G.E. Image 100 in this report) is a user-oriented 
system designed: for deriving useful information from remotely- acquired 
data, whether from satellites or aircraft. The G.E. Image 100 operates 
on the principle that all objects possess unique spectral characteristics or 
signatures. This signature uniqueness is used to "identify" similar 
features in an image by simultaneous analysis of signatures in multiple 
frequencies, such as the Landsat MSS. The G.E. Image 100 then assigns 
psuedo-colors or themes to those features having identical spectral 
characteristics. These themes may then be displayed on a color CRT, 
or super- imposed on the input image, or recorded. 
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Analysis of Landsat MSS Data 

Landsat MSS CCT's for snoivraelt periods during 1975 and 1976 were 
obtained for the portion of the Sierra Nevada that includes the American 

■| 

River basin. Snowpack extent during those two years differs widely. In 

1975, the snowpack was normal to slightly above normal. 1976 was the 

third-driest winter since 1900. MSS CCT’s were available for the dates 

given in Table 3. Only information from MSS bands 5 and 7 (0.6-0.7iJm 

and 0.8-l.lyin, respectively) were utilized in the G.E. Image 100 analyses 

because of the atmospheric haze problem in MSS band 4 (0.5-0. 6ym) and the 

of ■ ■ 

preference of the longer near-IR wavelength/MSS band 7 in place of the 
narrow, often saturated, band 6 (0.7-0.8ym) . 

The analyses entailed reading the data from MSS band 5 and 7 into 
the G.E. Image 100 with a resultant color display of the American River 
basin in the composited bands on a cathode ray tube monitor contained on tlie 
Image Analyzer Console. Through interaction with the Graphics Display 
Terminal, a histogram was generated summing MSS band 5 and 7 response, pixel 
by pixel, in a bivariate graph, MSS band 5 as the ordinate and MSS band 
7 as the abscissa. The pixels containing snow were grouped into three 
classifications as shown in Table 4. |The ranges for each classifi- 
cation were selected by viewing the CRT display for various groupings of the 
bivariate MSS 5 and 7 pixel distribution and then "alarming” the Image 
Analyzer Console television monitor to note which pixels had been 
designated a given category. In general, the reflectance limits for 
both MSS bands varied only slightly from scene to scene for a particular 
classification. In retrospect, a constant set of limits could have been 
specified for the classification of all scenes. If this had been done. 


TABLE 3. 

Dates 

River 

of available Landsat CCT’s covering American 
basin, 1975 and 1976. 

1975 


1976 

8 May 


14 April 

17 May 


23 April 

13 June 


20 May 

22 June 


29 May 

1 July 








TABLE 4 

RANGE OF REFLECTANCE IN ARBITRARY UNITS 
FOR VARIOUS CLASSES OF SNOW 





certainly the number of pixels in a given classification would change, 
but not sufficiently to alter the proportion of snow type to total 
snowcover. Furthermore, fixed limits would not be affected by changes in 
class limits that may result from changes in sun elevation. 

Images for the nine dates, with the three snow classifications are 

presented as Figs. Ha through Hi. These images were produced by a ■ 

Dicomed processor and show high elevation snow as yellow, middle elevation 

low 

snow as purple, and the/elevation snow as cyan. These three classifications 
are thought to represent fresh or nonmelting snow, melting snow, and partial 
snow cover, respectively. The remaining area contained within the American 
River basin, where vegetated, is reddish-brown; where bare rocks occur, 
the theme is gray. Gray pixels also seem to include some of the low 
elevation snow A comparison of the 8 May 1975 and 1 July 1975 images 
(Figs, lia and He) reveals that large areas changed from gray to 
red-brown, i.e., from snow to vegetation. By July 1 the largest areas 
of gray are limited to the high elevations in the vicinity of the Sierra 
crest . 

The percent of each snow type relative to the entire snow-covered 
area is provided in Table 5. Table 5 a^so compares computer- • 

generated SCA to SCA obtained using conventional snowraapping procedures 
practiced in the NESS (Schneider, et al., 1975). 

Comparison of the percentages of the three snow types reveals the 
lack of any consistent trend as snowmelt progresses. The percent of 
dry (yellow) snow decreases seasonally but erratically, with several 
reversals to larger percentages. Patterns for the categories of melting 
(purple) and partial (blue) snowcover are similarly erratic. These 
trends become more consistent, however, when the 
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Computer-generated snow classification for 17 May 1975. 




Fig. lib. 





Fig. 11 c. 


Computer-generated snow classification for 15 June 1975. 
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Fig. llh. Computer-generated snow classification for 20 May 1976. 
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Table 5, Comparison of Conventional Analysis vs. • Computer-Generated Analysis and Breakdown of Snowcover by Class as Pc 


Percent Snowcover Percent Snowcover Ratio 

Conventional Analysis Computed by G.E. Image 100 GB XdO/Conventional 


S May 1975 

45 

17 May 19 7S 

40 

13 Jun 1975 / 

16 

22 Jun 1975 

12 

1 Jul 1975 

7 

14 Apr 1976 

41 

23 Apr 1976 

26 

20 May 1976 

7 

29 May T976 

6 


FCHiDOUT ERAMi | 


Snowcover Breakdown by Type i| 
Yellow . Purple | 


0.60 

15 

1 

0.50 • 

7 

3 

0.31 

2 

<1 

0.17 

<1 

<1 

0.14 

. <1 

<1 

0.51 

13 

<1 

0.31 

1 

1 

0.14 

<1 

<1 

0.08 

<1 

<1 











by satellite, suggesting that sensor differences may exist between the 
MSS on Landsat-1 and the MSS on Landsat-2. Such differences would affect 
reflectance limits for a given snow classification and, possibly, to 
significant changes in the number of pixels placed in a given snow 
classification. 

Another problem is the omission of the contribution that gray pixels 
m^ke to the total snowcover Early in the snowmelt season when snow- 
free bare rock is at a minimum, most gray pixels appear to represent 
low elevation snow and should be included in the total pixel count for 
snow. However, inclusion of gray pixels as representative of snow in the 
late melt season would erroneously inflate the SCA, since most of these 
pixels would be representing bare rock rather than snow. It appears 
that the combination of MSS bands 5 and 7 is insufficient to differentiate 
between low-elevation snow and bare r 9 ck--at least for the American 
River basin. Inclusion of the percent of gray pixels in the American 
River basin would not solve the dilemma of fluctuation noted in Table 
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Surface Data 

A. Precipitation 

During the 1975 snowmelt period of May and June only light amounts 
of snow fell. The heaviest snovi/fall, 13 cm, occurred on 20 May at 
Twin Lakes (2386 m) in the southeast corner of the American River basin. 
All snowfalls occurred between Landsat overpasses and were melted 
before the subsequent satellite view of the basin. Non-snow precipita- 
tion (rain) was limited to the first week of May and to 24 June with 
generally less than 0.5 cm of rain occurring sporadically. Thus the 
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timing of precipitation events neither affected the retreat of the snow- 
line nor increased the surface reflectivity of the snow when examined by 
the MSS. 

Snowfall occurred on 14 April 1976 the same day as one of the 
Landsat overpasses. There is some cloud contamination in the southeast 
portion of the image (Figure lif), but a check of the NOAA-4 VHRR IR 1 
image (taken at approximately the same time as Landsat) failed to provide 
any information as to location or extent of the clouds. Evidently the 
clouds were low-level cumulus producing scattered, but at times, intense 
snow showers. One isolated cloud is readily apparent in the south-central 
portion of Figure Ilf ^ having a similar theme as the snow. Another cloud 
can be found within the snowpack, located near the east-central edge of the 
basin boundary. Though the cloud itself is not separable from the 
surrounding snow, the shadow stands out as an island of gray and blue in 
a sea of yellow. Other clouds if present over the snoivpack must be 
small, not casting discernable shadows. The overall cloud effect on the 
data displayed in Table 5 is thought to be small. | 

Except for the 14 April 1976 case, no occasional snowfalls are 
known. No precipitation effects in snow surface reflectance are noted 
for the remaining dates. 

B. Temperature 

Only eight stations in or adjacent to the American River basin 
report maximum and minimum temperatures to the National Climatic Center 
of EDS. One of these stations. Blue Canyon WSO airport, acquires 
temperatures every hour. Unfortunately Blue Canyon is not representa- 
tive of stations in the high Sierra, having minimum temperatures 5-15°C 


higher than other American River basin locations. Maximum temperatures 
are quite similar for all stations, including Blue Canyon, Elevations 
of the stations ranged from 1572 m (5156 ft) to 2386 m (7829 ft). 

For Blue Canyon, 1800 GMT temperatures (within 15 minutes of the 
Landsat overpass) average only 2,5°C less than the daily maximum, 
ranging from 0.5 to 4.0°C below the maximum ambient temperature. It is 
not known whether these 1800 GMT temperatures at Blue Canyon are 
representative of the temperatures at other areas in the American River 
basin for the same time. Copies of monthly thermographs for Twin Lakes, 
California, obtained from the Pacific Gas and Electric Company showed 
the 1800 GMT temperatures for the nine dates of Landsat coverage had 

climbed to 88 percent of the total diurnal rise in temperature. TemperaturW 
at Blue Canyon averaged 78 percent of the diurnal rise by 1800 GMT 
(1000 local) --perhaps a lower figure than all other stations in the 
American River basin. Thus an assumption of 80 percent of the diurnal 
rise by 1800 GMT for Sierra stations appears conservative and acceptable 
for this study. These estimated 1800 GMT temperatures (i.e. , minimum + 

0.80 (maximum-minimum) ) are presented in Table 6. No temperatures for | 
any of the nine dates are known to be subfreezing during acquisition of 
the Landsat data. In fact, only on 14 April 1976 are temperatures 
consistently estimated to be les? than 5°C at 1800 GMT. 

Ancillary temperature data were obtained from Dr. James Smith 
of the USDA Forest Service and Mr. Barney Clorety of the Sacramento 
Municipal Utility District (Table 6), Similarly these data indicate | 

subfreezing temperatures at 1800 GMT only on 14 April 1976 and then only 
above 2316 m. Upper air temperatures for 850 mb (1463 m) and 700 mb 
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TABLE 6 




STATION 


Esti mated S urface and Upper Air Temperatures at 

1975 

ELEVATION 8 May 17 May ' 13 June 
Meters Feet 


Blue Canyon WSO 

Donner Mem. St. Pk. 

Lake Spaulding 

Squaw Valley 

Tahoe City 

Tahoe Valley FAA AP 

Truckee R.S. 

Twin Lakes 
SMUD 1** 

SMUD 2 

SMUD 3 

SMUD 4 

SMUD 5 

SMUD 6 

USDA CSSL 


1609 5280 11.7* 21.7* 27.8* 

1810 5937 11.8 15.6 23.1 

1572 5156 12.7 19.9 26.8 

1900 6235 10.7 14,4 23.4 

1899 6230 9.3 12.8 20.9 

1906 6254 9.6 16.0 21.1 

1827 5995 

2386 7829 8.6 19.1 19.6 

1554 5100 18.6 25.2 31.2 

1798 5900 16.2 21.9 30.2 

2042 6700 13.0 19.8 25.7 

2316 7600 10.6 14.9 20.6 

2408 7900 14.9 - 26.5 

2682 8800 - 18.8 22.6 

2103 6900 12.7* 16.6* 21.7* 


850 mb 1463 4800 13.3 16.0 

700 mb ~ 3002 9850 1.0 6. 0 

* Actual 1800 GMT Temperature, Not Estimated. 
**SMUD: Sacramento Municipal Unility District. 



(3002 m) yield a similar finding (Table 6) . Thus it is likely that even 
at the Sierra crest (above 3050 m) the snow was melting during each 
Landsat overpass with the single possible exception of 14 April 1976. 
Comparison of Areal Snow Extent; Conventional vs. G.E. Image 100 Zoom 
Areal snow extent as measured by Zoom Transfer Scope (ZTS) or by 
computer via the G.E. Image 100, decreases seasonally as expected. However, 
the decrease in snow extent is more precipitous in the computer-derived 
values. The ratio of SCA via the G.E. Image 100 to the SCA via the conven- 
tional NESS methods varies from 0.6 for 8 May 1975 to 0.08 for 29 May 1976. 

A close inspection of the Dicoraed photographs provides at least one 
reason for part of the apparent rapid decrease in snowcover. Early 
during the spring melt (17 May 1975) the snowcover is more continuous, 
particularly at the high elevations (yellow) but becomes noticeably 
less continuous in the lowest elevation (cyan) . As the area of "nonraelting" 
snow decreases, the area of "discontinuous” snow grows. Thus near the 
end of the snowmelt season, only few patches of "continuous" snow remain. 

The G.E. Image 100 system tallies snowcover pixel by pixel, and is not 
concerned with the continuous snowline. The experienced snowmapper draws 
an enveloping, smooth continuous line, containing the SCA, but does not 
analyzing for the the individual discontinuities. Two conventional snow- 
cover analyses are shown in Figures 12a and 12b, and were completed using 
the Dicomed photographs for $olid lines) 17 May and 22 June 1975. Conven- 
tional snow-cover maps using standard 1:1,000,000 MSS band-5 images are 
shown in Figures 13a (17 May 1975) and 13b (22 June 1975) and have been 
added to Figures 12a and 12b (dashed lines) to facilitate comparison. SCA 
for 17 May 1975 was 40 percent and 33 percent, respectively, for conventional 
vs. Dicoraed (Fig. 13 vs. 12a) and for 22 June 1975 were 12 percent 
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Figure 12b. 


Comparison of SCA as determined from Fig. lid (solid line) 


with conventionally-derived SCA shown in Fig. 13b for 


22 June 1975 




and 9 percent, respectively, for conventional vs. Dicoraed (Fig. 13b) 
vs. 12b). All SCA values are well above the amounts achieved by G.E. 

Image 100 computer analyses, viz., 20 percent for 17 May 1975 and 2 
percent for 22 June 1975. 

Similar results were obtained for individual snow classifications. 

Two particular cases, 17 May 1975 and 22 June 1976, show conventional 
analysis of this "nonmelting, continuous" SCA (yellow) to be 12 and 2 
percent, respectively; values for the same dates and snow classification 
via G.E. Image 100 analysis were 7 and less than 1 percent, respectively — 
again significantly less than the conventional operational procedure. 
Discussion of Results 

A. G.E. Image 100 vs. Conventional Procedure as a Snow Mapping Tool 

Previous studies (e.g., Kirdar, et al.) have shoivn that the conven- 
tionally derived NOAA NESS SCA estimates agree well with ground observations 
Though SCA measurements were instantaneously achieved by electronically 
counting the classified snow pixels, the low values of SCA indicate the 
apparent difficulty of the G.E. Image 100 to properly define the various 
themes representing snow-covered area as reviewed by the MSS of Landsat. 

Most likely the discontinuous snow (blue) ia located in a mixture of 
coniferous forests and open areas, common throughout the Sierra Nevadas. 

The snow quickly melts from the trees, but often remains as deep as 2 
meters on the ground Unfortunately the forest partially obscures the 
Landsat MSS yiew of the snow-covered ground beneath. 

Another possibility for the observed differences in the tvtfo 

procedures used to calculate SCA could occur from the greater detail in 

area measurement obtained through the raultispectral analysis. There are 

snow-free areas within the snowpack which are too ntunerous and too minute 

to measure accurately with a compensating planimeter, but these small 
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melt areas or wind blotm areas can be detected and tallied with the 
automatic planimeter function of the G.E. Image 100. 

IVhatever the cause of the observed differences in the two techniques 
to map and measure snow, snow mapping is apparently best viewed synopti- 
cally using 1:1,000,000 MSS 5 images rather than the pixel by pixel 
(80-meter resolution) procedure follov;fd by the G.E, Image 100. The 
NESS operational snow mapping program described by Schneider, et al. 

(1975) has proved its usefulness as attested in communication received 
by real-time users of NESS-derived SCA values. For the present, manually 
prepared snow maps and resultant SCA values are timely and cost effective. 
Digital snow mapping using GOES data in near-real time is being done 
for six basins located in or near the Sierra Nevada, including the 
Anierican River basin (Tarpley, et al., in press). Plans to expand | 

the procedure to the remaining NESS basins are under consideration. 

Even if the G.E. Image 100 system produced equivalent snow maps matching 
those done by the conventional procedure, the technique would be more 
time-consuming (one must wait for CCT's) and costlier. 
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would be more time-consuming and costlier. 

B. Snow Classification using the G.E. Image 100 

Surface air temperature, interpolated to the time of Landsat over- 
pass, and upper air analyses indicate that the entire snovvpack was 
melting at 1800 GMT for all dates, with the possible exception of 
14 April 1976. On this date the freezing level may have been slightly 
lower than the 3050 m crest of the Sierras. The snowpack differences 
shown in the Dicoraed (colored) images cannot then be attributed to 
melting or non-melting snow. Rather the colors are more likely represen- 
tative of the number of hours of above-freezing temperatures the snow 
has encountered or the frequency of freeze-thaw cycles. Both of these 
parameters directly effect the grain size of the snow and thus the 
reflectivity of the snowpack (O'Brien and Munis, |1975). 

A second factor affecting measured MSS reflectances involves land 
cover classification. The predominent land cover at elevations where 
snow exists ranges from Ponderosa pine at the lower elevations, and 
fir and spruce in middle elevations, to Lodge-pole pine at the highest 
elevations. In the middle and upper elevations near and along the Sierra 
crest, large areas of sparse vegetation and bare rock predominate. In 
this region the Image 100 performs adequately, distinguishing snow from 
rocks when the snowline lies above 1800 m (5905 ft). At lower elevations 
early in the snowmelt season (see Figures lla-c; llf-g), the | 

classification scheme used omits the lowest elevation snoitfcover by 
categorizing the snow as bare rock. The effect of snow on Ponderosa 
pine, fir^ and spruce evidently results^ at Landsat resolution, in a land 
use pattern that appears so similar to bare rock that the combined 
information in MSS bands 5 and 7 is insufficient to differentiate between 
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snow and bare rock . Information contained in MSS bands 4 and 6 would not 
aid in the separation of snow from bare rock and thus the problem must 
remain unresolved using the present MSS system. 

. Clouds appear similar to snow (see Fig. Ilf) probably because ! 
of the low saturation level in band 5. Cloud shadows over the snowpack 
are classified as bare rock (see Fig. Ilf, 14 April. 1976) . . Snow 1 
classification images from the drought year of 1976 wrongly identify 
the bare, dried peripheries of the reservoirs exposed by low water 
levels as bare rock. 

An interesting observation is the limited amount of snow classified 


as "melting snow" (i.e., purple) in figures lla-i. Most of the snow classi-! 
fied as "partial cover" (i.e., cyan) should probably be placed in the 
"melting" snow category as well, especially as most of the partial- 
coVer snow lies at lower elevations than the "melting" snow (pusrple) . The 
8 May 1975 image (Figure 11a) is a good example- Ice House Reservoir ! 
at elevation 1657 m (5437 ft), located southeast of the largest ice- 
free reservoir (Union Valley Reservoir) in the American River basin, 
appears purple or classified as "melting" snow in its east and west 
ends. For a pixel to be so classified the snow must saturate MSS band 5 
and yet be a much lower reflector in MSS band 7 (reflectance not to 
exceed 44 on an arbitrary scale of 0 to 63) . Areas that are cyan 
maintain an identical range in MSS band 7 but are permitted to have MSS 
band 5 values as low as 70 on a scale of 1 to 127. Based on the range 
of MSS band 5 and 7 values for each classification and knowing that 
temperatures were above freezing everywhere in the basin at the time of 
the satellite overpass, it appears to be more appropriate to classify 
purple snow as middle metamorphosed snov< and cyan snow as late 
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meitamorphosed snow (or initially melting and actively melting snow, 
respectively). Yellow snow would be better labeled early metamorphosed snow 
rather than fresh snow since the date of last snowfall may have occurred 
several weeks prior to the satellite overpass. These classifications are 
purely arbitrary and without sufficient ground truth cannot be substan- 
tiated. However, these classifications appear to be more realistic than 
the ones proposed earlier in this section. 

The apparent reason for the small number of pixels colored purple is 
the result of only one value being acceptable for MSS band 5, viz., 127 or 
saturation. This restriction permits only 31 possible combinations of MSS 
band 5 and 7 reflectances for the purple snow category. Possible combi- 
nations for yellow and cyan snow are 780 and 1,922, respectively, or 25 
and 62 times as many opportunities to have a yellow or blue pixel than 
one that is purple. Table 7 presents the ratios of snow classified | 

as "Yellow" and "cyan" to the snow placed in the "purple" category. Only 
once, 14 April 1976, when the lowest surface temperature of all nine cases 
occurred, did the measured ratios exceed those expected by chance. In 
other words, there were generally more purple-alarmed pixels than xvould 
be expected to occur in a purely random process. 

Similarly, when comparing yellow and cyan classifications, the 

occurrence of yellow pixels in proportion to cyan pixels exceeds that 

expected by chance, except for the 23 April 1976 case. Thus the cyan 

snow classification is less common in comparison to chance and to both 

the purple and yellow snow classification. Statistics are easily 

generated; however, to grasp the real significance of the thematic 

Approach to snow separation by type is much more difficult. Without 
Sufficient ground truth and because of the saturation of the MSS by 
snow, the thematic approach using Landsat MSS data generally cannot 
provide meaningful results. 
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■TABLE 7. 


^ Snow TyRe__and--Com ,oJ_Yellow. iind Cyan-Sno>.v-to Purple Snow in Percent 


8 MAY 1975 17 MAY 1975 15 JUNE 1975 22 JUNE 1975 1 JULY 1975 

Ratio Ratio Ratio Ratio Ratio 

No. of Snow type: No. of Snow type; No. of Snow type: No. of Snow type: No. of Snow type: 

Pixels purple snow Pixels piorple snow Pixels purple snow Pixels purple snow Pixels purple snow 


17 MAY 1975 
Ratio 


22 JUNE 1975 
Ratio 


1 JULY 1975 

Ratio 


Purple . 10,389 
Yellow 156,852 


15.10 


29,411 

80,022 


18,544 


1.84 


Cyan 116,197 


11.18 108,036 





14 APRIL 1976 23 APRIL 1976 

~~ ! Ratio Ratio 

Snow No. of Snow type: No. of Snow type: 

Type Pixels purple snow Pixels purple snow 


> 

Purple 

292 

1 

12,208 

1 

647 

1 

877 

1 

in 

Yellow 

155,079 

462.60 

15,833 

1.30 

5058 

7.82 

1071 

1.22 


Cyan 

91,148 

312.5 

60,624 

4.97 

6005 

9.28 

3299 

3.76 



20 MAY 1976 
Ratio 

No. of Snow type: 
Pixels purple snow 


29 MAY 1976 
Ratio 

No. of Snow type: 
Pixels purple snow 



\ 

I 

r 

i - 


r 





Summary 


From G.E, Image 100 computer tapes, Dicoraed color images of 
Landsat scenes were generated. Various colors, each assigned to a 
bivariate range of MSS band 5 and 7 arbitrary reflectances, were used 
thematically in an attempt to classify the changing condition of the 
American River basin snowpack. From the histogram of MSS 5 and 7 values 
three classifications or themes vvere readily apparent and were thought 
to repi'esent partial snowcover (cyan), melting snow (purple), and fresh 
or nonmelting snow (yellow). Inspection of available surface and upper 
air temperatures revealed that melting conditions existed on all Landsat 
overpasses with the possible exception of 4 April 1976. On this date, 
it is likely that the highest elevation snow surface (above 2500 m) was 
at or slightly below freezing at the 1800 GMT satellite acquisition. 

Rather than separate melting from non-melting snow, the thematic 
approach has apparently divided the snow into various degrees of metamorphisra 
Possible classifications of these three categories could be actively 
melting or late metamosphosed snow as cyan, early melting or middle 
metamorphosed snow as purple, and initially melting snow or early 
metamorphosed snow as yellow. Unfortunately snoivpack ground truth 
information in this basin was sparse to nonexistent and precluded any 
precise evaluation of the snowpack condition. 

Another shortcoming noted in the G.E. Image 100 thematic separation 
was the inability to note a spectral difference between lowest elevation 
snow and bare rock. This inability is most critical late in the snoi<- 
melt season when the snowline has receded sufficiently to expose large 
areas of bare rock. At this time of year, classifying bare rock as low- 
elevation snow would lead to inflated values for total snowcover. 
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The G.E. Image 100 computer system is quite capable of handling 

■ ^ ! - ■ ' ■ 

the sophisticated multispectral techniques necessary to develop themes 
for snowpack classification. Unfortunately the Landsat MSS sensor 
limits the G.E. Image 100 's usefulness because of low energy-level 
saturation and inadequate multispectral information. Not until the 
MSS system is improved by eliminating these problems will substantial 
progress result in a worthwhile assessment of snowpack parameters. 

Saturation limits of future HSS's must be raised to avoid the 
currently-occurring saturation of the MSS when imaging snowfields. 

The AVHRR on TIROS N, NOAA's new operational polar-orbiting satellite, 
has saturation limits that exceed the reflectance observed from snowpacks 
and demonstrates thas sensor saturation need not be a problem. 
Preliminary research has been completed regarding the spectral response 
of various natural snow samples in the visible and near-IR (O’Brien and 
Munis, 1975). Research at CRREL is currently being conducted, in situ, 
on the effects of various substracts on snow reflectance. These 
new data combined with CRREL 's previous study should provide insight 
into which visible and near-IR spectral regions provide useful 
information concerning snowpack properties. 
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1978: Snowfield Assessment 
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Snow in Forest Areas 


Introduction 

Snow mapping in dense forested areas, whether deciduous or coni- 
ferous, was recognized as difficult from the inception of our research 
into snow mapping at NESS. A summary of the findings of a report pre- 
pared for NESS by Eschner et al . (1977) investigating the problem of 
snow monitoring in forested areas via remote sensing follows. 

Ob j ectives . 

The study (Eschner, et al., 1977) comprised three principal 
efforts: (1) a quantitative consideration of masking effects resulting 

from Vegetative and topographic characteristics, (2) a detailed analysis 
of snow reflectance under vegetative cover using Landsat MSS and NOAA 
VHRR data, and (3) experimental mapping of snow extent via manual, 
digital, and electronic methods. Data used included NOAA digital VHRR 
visible and thermal tapes, Landsat digital MSS tapes, ground observations 
of reflectance, and low-altitude photography. 

Data 

The handling of the NOAA digital VHRR data required the development 
of various software programs; a description of the five required programs 
is provided in the report (Eschner, et al., 1977, pp. 18-21). Based on 
techniques described by Algazi, et al. (1976) and Legeckis and Pritchard 
(1976), the VHRR data were corrected for geometric and panoramic 
distortions as well as the problem of skew. Programs for processing 
digital Landsat data were developed concurrently with the NOAA software. 
Data Analysis 

A correlation analysis between NOAA VHRR visible brightness and 
site characteristics such as percent conifer, slope and roughness, 
slope and shading, and snow depth, as tabulated in Table 8, showed 
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Table 8 


Correlation Analysis 
al., 1977) 



Correlation Coefficients, r's 

Snow Depth 
Percent Conifer 
Slope/Roughness 
Slope/Shading 





s 


over 

Snow Survey Sites 

(After 

Eschner, et 

NOAA Visible Brightness, expressed as: 
Raw Corrected % tile Norraal'd 

.15 

.06 

.03 

.10 

o 

CM 

-.29 

.27 

.17 

.22 

- . 14 

.23 

.21 

• 

o 

00 

.01 

.07 

.07 
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only weak correlations. Such weak correlations did not warrant the use of 
a regression equation to estimate snow depth based on site characteristics 
and visible brightness. The lack of correlation is apparently caused 
by the poor signal -to-noise nature of the VHRR data. 

Table 9 statistics demonstrate that the smaller range of the ! 

Landsat MSS response suggests greater sensitivity and better signal-to- 
noise ratio than the NOAA VHRR visible band. 

Contrary to earlier findings of Wiesnet et al. (1975)^ who reported 
a decrease in reflectance of snow in the Cascade Mountains as a function 
of increasing vtfavelength, Eschner, et al. (1977) in the Adirondack Mountains 
measured a maximum response (i.e., reflectance) in the near-IR. Table 
10 presents a summary of the grouped MSS data (for deciduous, conifer, . ‘ 
and open areas) , first as raw data (digital counts) then in terms of 
radiance (pW/cm sr) . The findings are graphed in Figure A. 4.1 and show 
in general an increase in radiance for all types of ground cover with 
increasing wavelength. The only exception is the result for deciduous forests 
on 25 February 1976 where the band 6 radiance of 2.18yW/cm^sr fell below j 

2 

the band 5 value of 2.32pW/cm sr. Eschner et al. (1977) attribute the 
increased reflectance of forested areas "to the high density of vegeta- 
tion present in the study area" (adjacent to Cranberry Lake) and in the 
northeast in general. This explanation may have credence over 
coniferous sites, but would not apply to deciduous forests where no apparent 
explanation exists for the noted increased radiance. 

Hard-copy visible-band NOAA VHRR imagery, analyzed following the 
procedure described by Schneider et al . (1976) , yielded snow extent 
maps in close agreement with ground surveys. Some difficulty was noted 


Table 9 


Satellite Imagery Data Ranges, 1976 Dimensionless Values, 

Absolute Scale of 0 to 255. (After Eschner, et al., 1977). 

96 Percentile Look Angle 




Limits 

Range 

Mode 

Mean 

St. Dev. 

a 

NOAA 

15 Jan 

183-224 

41 

205 

204 

11 

.328 

(visible) 

17 Jan 

191-233 

42 

223 

220 

11 

,280 


19 Jan 

192-231 

39 

223 

224 

12 

.210 


23 J an 

201-230 

29 

215 

225 

11 

.096 


23 Feb 

196-231 

35 

223 

228 

14 

.144 


25 Feb 

186-228 

42 

215 

227 

18 

.203 


24 Mar 

188-226 

38 

215 

229 

19 

.371 


7 Apr 

209-230 

21 

223 

220 

6 

.057 


9 Apr 

209-232 

23 

223 

220 

6 

.112 

Lands at 

25 Feb 

23-177 

154 

43 

66 

39 


(band 5) 

24 Mar 

20-153 

133 

45 

60 

31 


All units 

are digital 

counts on 

a scale 

of 0-255. 
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and Radiance Estimates. (After Eschner, 

t Radiance Estimate 

(yw/cra^sr) 

Band 

7 4 5 6 7 


61 

2.60 

2.32 

2.18 

3.87 

329 

40 

1.-57 

1.21 

1.37 

2.54 

95 

82 

1.43 

1,33 

0.99 

1.29 

315 


54 

2.28 

1.78 

2.01 

3.43 

298 

41 

1.22 

0.76 

1.29 

2.60 

168 


# pixels 
sampled 


40 1.15 0.87 0.70 0.63 


240 


INTENSITY. DIGITAL COUNT 




(0 

I 1.0 



4 5 6 7 

MSS BAND 


2 0.6 

< 

O 

< 

OC 0.4 

O 

UJ 

is 0.2 

s 



LEGEND 

DATE 

2/25 3/24 


4 5 6 7 

MSS BAND 


O CONIFER 
□ DECIDUOUS 


Graph showing Landsat brightness levels and radiation 
estimates o£ snowcover in deciduous and conifer areas 


of the Adirondack Mountains, (after Eschner, et al.. 



in the analysis o£ 1976 images prior to late March when the snow-covered 
Adirondack region consistently appeared as dark as the non-snow areas 
to the south. 

Digital analysis of the VHRR data from 7 April 1976 required an 
averaging technique to remove noise in the data. Once the noise- induced 
"salt and pepper" tone was removed, the resulting product resembled 
snow maps produced using standard VHRR imagery and the Zoom Transfer 
Scope. The digital analysis technique for snow mapping is more 
advantageous than conventional hard-copy photograpic interpretation in 
at least two aspects. The imagery may be geometrically corrected for 
precisely relating image data to watershed boundaries, and groiind 
observations may be more accurately located on an image to assist in 
determining the threshold level separating snow from non-snow areas. 
Summary 

NOAA visible VHRR digital data present little information on snow 
conditions in the Adirondacks during winter and early spring. This 
lack of information appears to be caused by a poor signal-to-noise 
ratio in the data. After mid-March, discrete areas of high reflectance 
correspond well to areas of snowcover indicated by ground surveys, 
providing timely information to the hydrologist during the spring run- 
off period. Digital generation of snow maps from VHRR visible data is 
hanipered by noise in the data, but the application of a simple averaging 
technique was found to greatly increase the utility of the data. This 
procedure is more precise, but more costly than either the manual or 
computer-assisted interpretation of snow from hard-copy NOAA imagery. 

The increased near-infrared response observed in Landsat MSS 6 
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and 7 from dense coniferous and deciduous forests warrants further 
examination. Future spectral studies involving the discrimination of 
fresh and metamorphosed snow should consider the effects of varying 
vegetation type and density as a function of resolution element size. 


References ; (Section on Snow in Forest Areas) 

Algazi, R., 1976: An Integrated Study of the State of California using 

Remote Sensing Techniques, Annual progress report. University 
of California, Davis, pp. 58-70. 

Eschner, A.R., Lillesand, J.M., and Meisner, D.E., 1977: Satellite 

Remote Sensing of Snowcover in the Adirondack Mountains, Final 
Report, 88 pp., NOAA/NESS Grant No. 04-5-158-43. 

Legeckis, R. and Pritchard, J., 1976: Algorithm for Correcting the 

VHRR Imagery for Geometric Distortions due to the Earth's Curva- 
ture and Rotation, NOAA Tech. Memo. NESS 77, 15 pp. 

Schneider, S.R., McMillan, M.C., and Wiesnet, D.R., 1976: River Basin 

Snow Mapping at the National Environmental Satellite Service, NESS 
Tech. Memo. NESS 83, 27 pp. 

Wiesnet, D.R., McGinnis, D.F., and McMillan, M.C., 1975: Evaluation of 

ERTS Data for Certain Hydrological Uses, Final Report to NASA/GSFC 
Contract No. 432-641-14-04-03, 92 pp. 


A-64 


Snowcover Runoff Relation in American River Subbasins 

I. Introduction 

The monitoring of changes in the extent of snowcover began with the 
advent of operational satellites. The feasibility of using snowcover 
changes to extract useful information regarding resultant snowmelt 
runoff was reported by Range and Salomonson (1977), Saloraonson and 
MacLeod (1972)» and the U.S. Department of Interios (1974). 
of Range and Salomonson involved correlating the snowT covered area (SCA) 
at its greatest extent, viz., April 1-20, (just before the beginning 
of the melt season) with seasonal runoff volume for April-July. Two 
watersheds in the Himalayan region were studied, yielding correlations 
of 0.91 and 0.94 (r) for straight line best fits. These high correla- 
tions must be tempered by the fact that only six and eight years of 
data were available for the two regressions. i 

II. Description of Subbasins 

A natural extension of the snowmelt-runoff studies is to shorten 
the period of runoff forecasts, for example, reducing 3-raonth volumes 
to monthly or perhaps weekly or daily runoff estimates based in part 
on changes in satellite-deserved SCA. Because we have gathered a 
considerable amount of Landsat data on the American River basin, it was 
decided to concentrate on subbasins of the American River. Only three 
subbasins of the American River have gaging stations on streams whose 
flows are unregulated upstream of the gage. These stations are Pilot 
Creek above Stumpy Meadows Lake (11431800), drainage area 30.3 km^ 

(11.7 mi^); Duncan Creek near French Meadows (11427700) , drainage 
area 25.74 kra^ (9.92 rai^); and North Fork of Middle Fork American 


River near Foresthill (11433260), 230.3 km^ (88.9 mi^). The higher 
elevations of each of these subbasins receive snow during the winter 
and spring seasons --though the snow depth and cover within each 
subbasin vary greatly from year to year. The range in elevation is 
1606 to 2195m for Duncan Creek, 1304 to 1725m for Pilot Creel^ and 396 
to 3189m for North Fork of Middle Fork American River. None of the 
gages predate 1960 and range from 13 to 18 years in length of record. 

Figure 15 shows the location of the three subbasins within the j 

American River. Figures 16a, 16b, 16c, show details of • 

Duncan Creek, Pilot Creek, and North Fork of Middle Fork American 
River basins » respectively. 

III. Data and Data Analysis 

Average daily runoff values, in cubic feet per second, are 
published yearly by the U.S. Geological Survey in cooperation with 
individual state agencies (USGS, 1973-76) . The 12-month period of 
record begins on the first of October and is known as the water year. ) Thus the 
1973 water year covers the period October, 197^ through September, 

1973. Plots of the daily discharge for the three subbasins for the 
year 1973, 1975, and 1976 are presented in Figure 17a for Duncan 
and Pilot Creeks, and in Figure 17b for North Fork of Middle Fork 
American River, respectively. Satellite-derived SCA in terms of 
percent of basin are listed in Table H« Note that Landsat data i 

are available for only a few dates per year, only *;wice for the 
drought year 1976, To compensate for the paucity of Landsat- derived 
SCA, NOAA VHRR SCA data were combined with the Landsat estimates thus 
providing a more complete record. SCA regressions were performed 
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Figure 15. jSketoh map showing the location of the Duncan Creeki 
Pilot Creek, and North Fork of Middle Fork American 


River subbasins. 



Figure Map of the Pilot Creek subbasin. 
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Daily discharge hydrogrs^phs of Pilot Creek and Duncan Creek 
for March -June 1973, 1975, and 1976. 
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Figure 17b. 


Daily discharge hydrographs of North Fork of Middle Fork 


American River for March-June 1973, 1975, and 1976. 




Table 11 • Snow Cover Area of American Rivex' Subbasins in Percent for 
1973, 1975, 1976 using Landsat MSS Band 5 Imageiy 
AMERICAN RIVER SUBBASINS 


Date 

N.F. of M.F. 
American River 

Duncan Creek 

Pilot Creek 

4/3/73 

48 

85 

64 

4/21/73 

40 

93 

48 

5/9/73 

19 

92 

24 

S/27/73 

4 

62 

0 

4/11/75 

74 

91 

84 

5/8/75 

40 

93 

73 

5/17/75 

25 

86 

22 

6/4/75 

4 

67 

0 

4/14/76 

29 

75 

69 

4/23/76 

8 

47 

0 
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relating NOAA-derived and Landsat -derived American River SCA as well as 
regressing individual subbasin SCA's against the American River SCA's 
(Landsat data only). Even though the data were sparse, meaningful 
relationships were found for Pilot Creek and the North Fork (NF) of 
the Middle Fork (MF) of the American River and are given below: 

NOAA vs. Landsat, American River Basin 


1973: 

y 

= 1.18 + 

l.Olx, 

r = 0.99 

y = Landsat-derived SCA 

1975: 

y 

= 6.99 + 

1.22x, 

r = 0.99 

X = NOAA-derived SCA 

1976: 

y 

= 0.81 + 

X 

00 

o 

r = 0.99 


N.F. of M. 

F. 

American 

River 




1973: y =-29.00+1.56 (1.18 + l.Olx) = -27.16 + 1.58x 

1975: y =-29.00 + 1.56 (-6.99 + 1.22x) = 39.90 + 1.90x 

1976 y = -29.00 + 1.56 (0.81 + 0.98x) = -27.74 + 1.53x 

X = NOAA-derived American 
River basin SCA, 

REPRODUCIBILITY OP Ti . in percent 

ORIGINAL PAGE IS POOR 

y = N.F. of M.F. American 
River basin SCA, 
in percent 

Pilot Creek 

1973: y = 21.90 + 1.71 (1.18 + l.Olx) = -19.88 + 1.73x 

1975: y = 21.90 + 1.71 (-6.99 + 1.22x) = -33.85 + 2.09X 

1976: y = 21.90 + 1.71 (0.81 + 0.9Sx) = -20.51 + 1.68x 

X = NOAA-derived American. 
River basin SCA, in 
percent 

y = Pilot Creek SCA, in 
percent 

Data for Duncan Creek were too scattered and showed no apparent 
relationship to American River SCA. 
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A thorough but somewhat outdated study of the snowmelt process 
is presented in Snow Hydrology (1956). The study indicates that 
temperature is a good index to snowmelt, to be explicit, the number of 
degrees that the maximum temperature exceeds the freezing point. Leaf 
(1975) has incorporated SCA into a computerized short-term strearaflow 
forecast model, that also includes the residual snowpack water equiva- 
lent as an independent variable. For our study of the American River 
basin, several independent parameters were regressed in an attempt to 
estimate daily runoff flows (in cfs) . The independent variables 
include SCA (for the American River basin and for the individual sub- 
basin), the number of degress the daily maximum temperature exceeded 
freezing in degrees F, the snow on the ground (inches), snowfall 
(inches), the water equivalent of the snowpack (inches) , and the amount 
of precipitation (in inches of water) . A final independent parameter 
was Julian day. All data were available on a daily basis, except SCA. 
Missing values of SCA were estimated by linear interpolation when rather 
uniform maximum temperatures occurred and precipitation was absent. 
Adjustments in the estimated SCA were made to account for snowfalls, 
widely varying temperatures or large rainfalls. Most adjustments 
involved less than several percent change in SCA. Estimates of SCA 
were made taking into account as many meteorological variables as 
possible; nevertheless some subjectivity was introduced in the 
decision process. The daily parameters, including the dependent 
variable, runoff, were punched onto computer cards and run through 
a packaged stepwise, multivariate regression on the NESS IBM 360/195. 

Because of the high correlation between SCA for the American River 

. . • , ( . 

basin and the two subbasins (0.94 to 0,99), a second run was made 
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exchanging SCA values. Data were available for 1973, 1975, and 1976; 
but no regression analysis were performed for Pilot Creek for 1976 
because of the rapid melt of the snowpack and the very low, almost 
constant, runoff. 

IV. Discussion of Resu lts 

Multiple regression correlation coefficients C^'s) ranged from 
0.74 to 0.96. The regression analyses usually included only three of 
the seven independent variables; the others contributed little to 
reduce the variance (suras of squares). Of the 10 analyses performed, 
only two required more than three variables to reach the maximum 
reduction in the variance. For these cases, five variables were used. 

Two other analyses required only two variables to reach maximum r. 

Snow-covered area* by providing the largest reduction in 
the sums of squares, was the most important variable in six of 10 
runs, occurring as the second and third most important parameter twice, 
respectively. Precipitation and day of the year were each most 
efficient in reducing the suras of squares twice. Snow on the ground 
appeared as only the second-most useful variable in the analyses, but did 
so in eight of the ten runs. In third place, in addition to SCA, were 
day of year (four times) and precipitation, twice. For the two cases 
when five variables were used, snowfall occurred once in both the fourth 
positions,' water equivalent of the snowpack once in fourth 
position, and day of year once in fifth place; 

The occurrence of SCA as an important variable in all the 
regression analysis is noteworthy. Further, the SCA for the American 
River basin in four of the five cases reduced the variance slightly 
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more than did the SCA for the individual subbasins. The reason ’for 
this slightly greater reduction is not known, but is of little 
significance because of the high correlation between the two variables. 

It suggests that one need not rely on the higher resolution, less frequent 
Landsat imagery to provide SCA. Rather»the lower resolution NOAA data 
should suffice in estimating runoff of small watersheds (as small as 
30.3 km^) as long as the SCA of the subbasin correlates well with the 
larger, NOAA VIlRR-size river basin. 

The other useful variables are standard meteorological or 
hydrologic parameters, and though some are not acquired at many 
observation stations (such as snow on the ground and snowfall), these 
variables are usually available near, if not within, the basin perimeter. 
In this study. Blue Canyon Weather Service Office (WSO) supplied the 
snow parameters ndd is located approximately 75. km from the subbasins. 

There is an indication that the parameters used are not always 
well -correlated i\fith runoff. For both basins, 1975 regressions 
have very low r’s of 0.74 and 0.75. For this year the Pilot Greek 
SCA was only the third most important parameter for reducing the sum 
of squares. SCA i>ras still the most effective parameter for the NF 
of the MF American River, but accounting for an r of only 0.54. 

Obviously other factors affected the runoff during 1975--factors not 
included in the chosen parameters. A check of the hydrographs for 
the NF of MF American River show 1975 to be rather irregular with 
three distinct peaks during the snowmelt season beginning in late 
March. 1976 shows a much different hydrograph response having only 
one double-peaked maximum with a long, almost smooth, recession 


curve. It vvas 1976 that the correlation coefficient was the highest 
and was accomplished using only the day of year and SGA variables. 

V. Summary 

In summary, a simple regression model of readily available 
meteorological parameters, satellite-derived SCA, and day of year pro- 
vides estimates of daily average runoff with r's in excess of 0.89. 
Such estimates are obtained during snoivmelt seasons where the hydro- 
graph is very regular, i.e., having one well-defined peak followed 
by a rather smooth, recession curve. During snowmelt seasons 
punctuated with many precipitation events, the regression model 
performs poorly. Thus SGA appears to contain information for those 
who are interested in better quantifying daily runoff in the American 
River basin and possibly for quantifying runoff in other Sierra 
Nevada basins similar to the American River. The potential to 
predict runoff, indirectly, in much smaller subbasins, the size of 
30 km^, in the American River using SCA was demonstrated. 
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SOIL MOISTURE STUDIES 


Introduction 

The need for improved soil moisture estimates is demonstrated by 
the crude approximation currently Used for hydrologic design (USDA, SCS, 
194:7) and streamflow estimates (Burnash, et al, 1973) . Most soil 
moisture values are actually based on precipitation values from rain 
gages that are weighted by a decay function of time or on accumulated 
5-day antecedent rainfall. Soil moisture measurements are rarely taken 
in situ, and when taken individually are point values, undoubtedly 
unrepresentative of the entire river basin. There is a real need for 
areally averaged soil moisture values that reflect the condition of the 
watershed aud its response to rainfall more accurately. 

Basin-wide soil moisture estimates are best obtained from remote 
platforms such as aircraft or satellites. Sensors scanning in such 
portions of the electromagnetic spectrum as the near-IR (Merritt and 
Hall, 1973), thermal (Idso, et al., 1975, Jackson, et al*, 1976), gamma 
ray (Feimster, et al, 1975), and microwave (Schmugge, 1976; Schmugge, 
et al., 1976a; Schmugge, et aL, 1976b) have been used in an attempt to 
obtain basin-i^ide soil moisture. An extensive study of active microwave 
work with soil moisture has been done by Batlivals and Ulaby (1977), 
while work by Schmugge, et al. (1976a) covers the passive microwave area. 
Unfortunately the resolution of satellite-borne passive microl^ave data 
is currently limited to approximately 25 km. 

Field studies were conducted near Phoenix, Arizona, and near 
Luverne, Minnesota, as part of this Landsat-2 investigation. The 
objective of the investigation is to monitor soil moisture at NOAA test 
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sites and to compare ERTS (Landsat) spectral reflectance measurements 
with not only ground truth (soil samples) but with concurrent gamma-ray 
aircraft' data collected and furnished by ongoing in-house NOAA research, 
and with concurrent NASA PMIS overflights. 

Phoenix Test Site 

Description of Site and Overview of Experiment 

The test site in the Phoenix area was selected for several reasons: 


% 


' 

t ■ 


I 

I 

U 


1) it is irrigated land, thus affording a variety of soil moisture; I 

2) the fields are regular and set along straight lines oriented north- | 

■ ; , ' . I 

south; 3) it had been studied previously in 1974 (Schmugge, et al, 1976b); : 

and 4) cloud-free conditions are relatively frequent. The experiment was i 

timed to coincide with that time when a large number of bare fields were | 

were . ^ ■ . I 

either planted in cotton or / being prepared for planting; therefore . ^ 

soil moisture conditions would vary widely. Figure 1 shows the general ' 

. i 

location of the test sites, and figures 2a-c present detailed maps of | | 

each field and the location of soil samples in each field. 1 

The experiment involved participants from NASA/JSC and GSFC, USDA/ 

SCS, NOAA/NESS, EG§G and Bittinger and Associates, as well as students 

from Arizona State University. Various remote sensors collected data in the I'garama 

ray, microwave, visible, thermal and near-IR portions of the speetrum for 

correlation with the soil moisture values determined from soil samples 

taken in the fields. Table 1 summarizes the remotely sensed data 

sources. 

The goal of the combined experiment was to evaluate the usefulness 
of remote sensors for monitoring soil moisture. The major thrust centered on 
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TABLE 1. Sununary of remote sensors used to collect data of Phoenix 


Wavelength 

test site. 
Sensor 

Platform 

Supplier 

Date 

Collected 

Altitude 


0.5-l.lpra 

MSS 

Landsat 

NASA 

3/17/75 

900 km 


10. 5-12. 5pm 

VISSR 

SMS-2 

NOAA 

3/18/75 

37, 500 km 


2. 8 cm 

PMIS 

P-3 

NASA 

3/18/75 

o 

o 

• 

CO 

km 

21 cm 

MFMR 

P-3 

NASA 

3/18/75 

0.5 § 0.8 

km 

lA (10~^lcm) 

Nal(Tl) 
(gamma ray) 

Beechcraft 

EG8G 

3/18/75 

0.15 km 



TABLE 2. Net and ^®^t 1 gairana photopeak count rate from the 

March 18, 1975 survey. (From Feimster, et al., 1976) 


Survey Spectrum 

Time 

of 

Net Gamma Count Rate in Photopeak (cpm) 

Line* 

No. 

Day 

4Uk 


lA 

203 

A'l 

2397 

534 

IB 

207 

AM 

2777 

592 

1C 

211 

AM 

2785 

598 

(Average values 

for survey 

line 



segment) 



(2673) 

(577) 

lA 

225 

PM 

2734 

573 

IB 

22e' 

PM 

2621 

583 

(Average values 

for survey 

line 



segment) 



(2628) 

(585) 

(Average values 

for i;v’hole ^ 

line) 

(2651) 

(581) 

*The symbols A, 

B , and C refer tc 

) the successive 

passes over line 1. 


h 

h 
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microwaves (NASA) with secondary efforts involving gamma ray, near-IR 
and thermal (NOAA) data. 

Soil Sampling and Soil Moisture Data 

Sampling was done at four points in each 16.9-hectare field selected. 
Points were located approximately one-fourth the width of the field 
from the adjacent boundaries (see sampling locations, figures 2a-c) and 
are thus at or near the centers of each quadrant of the field. 

Primary tillage practice was listing, i.e., the furrowing of 
the soil in which row crops are soivn in the furrows. Most fields, bare 
particularly, had listed furrows where irrigation water is distributed 
along the furroiv, and moisture reaches the listed ridges by infiltration 
and capillary movement. At each data point five samples were collected 
on the tops of ridges and five were collected at the bottom of furrows. 

The ten samples collected at each location were as nearly representa- 
tive as physically possible of the 0-1 cm, 1-2 cm, 2-5 cm, 5-9 cm, and 
9-16 cm depth zones. These sampling depths were chosen by the NASA/ 

Texas A§M team to test theoretical models of microwave emission. Weighted 
composite averages were calculated to represent composite soil moisture 
from the surface to the bottom of each zone to be used in the development 
of en 5 )irical models for determination of soil moisture itfith microwave 
equipment . 

A report by Blanchard (1975) px'ovides details of the sampling tech- 
niques and sampling processing. Appendix A to his report gives all soil 
moisture values including various soil moisture averages for different 
depths and for entire fields . Location of the fields and soil temperature 
data are also provided. 
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Microwave Study 


The relation of passive microwave to changes in soil moisture has 
been studied at the Phoenix test site each year from 1972 to 1975 
(Schmugge, et al, 1976a; 1976b; and Schmugge, 1976). The results of the 
various experiments, culrainating in 1975, showed that radiometers 
operating at 21 cm have good sensitivity to soil moisture variations, 
about a 50°K change as soil dries from field capacity. The sensitivity 
is maintained even for a moderate vegetation canopy. For wavelengths 
less than 3 cm, vegetation effectively screens the soil surface thereby 
limiting the usefulness of microwave radiometers for soil moisture 
sensing. Further study is needed to determine whether the additional 
polarization obtained at off-nadir orientations is worth the sacrifice 
of poorer spatial resolution. For a complete review of the microwave 
investigations one should consult the above referenced reports. 

Gamma Ray Study 

The purpose of the survey was to determine soil moisture by detecting 
differences in the terrestrial gamma flux attenuation as a function of 
soil moisture content. Most of the gamma radiation from the naturally 
occurring radioelements in the soil, measured above ground, originated 
within a few centimeters of the surface. The presence of moisture in the soil 
causes an effective increase in soil density resulting in an increased 
attenuation of the gamma flux for wet soil and a corresponding lower flux 
above the ground surface. Selected samples of the surface soil, analyzed 
gravimetrically, provide the necessary calibration data to allow a 
quantitative determination of aerial measurements. The EG§G aerial 
measurements of terrestrial passive gamma radiation were performed 
simultaneously with the NASA microwave measurements and the ground 


sampling, and offers an independent set of data from which to establish 
correlations and evaluate the NASA microwave data. 

The survey lines along which soil moisture measurements were made 
(25 km) C375u) 

consisted of tvtfo 16-mile/ long lines situated 1,200 feet/ on either side 

of 91st Avenue in the western section of Phoenix as shown in figure 1. 

These lines were oriented in a north-south direction and covered flat 

farralcind, most of which was unvegetated. The character of the region 

surveyed is illustrated in more detail in figures 2a- c, with indications 

of the location of ground-based soil sampling points. 

Equipment 

The EG^G aircraft survey consisted of measuring the natural 
terrestrial gamma-ray emission using an unshielded fourteen crystal array 
of 4-in by 4-in Nal(Tl) detectors. The gamma radiation comes primarily; 
from the naturally occurring in the surface soil, together with 
^^^Bi and ^®^T1 from the uranium and thorium decay chains. The aerial 
radiation measurements included the collection of gross gamma counts and 
gamma energy spectra. The gross gamma-count system consisted of an 
amplifier discriminator computer unit that counted and recorded the 
total number of gamma rays of energy greater than 50k eV that were 
detected during a one-second interval. The gross gamma-count rate 
(number of gamma rays detected per second) v^as digitally recorded along 
with aircraft position and altitude every second. As a backup and 
complement to the digital recording of the gross -count data, a record 
was made on a continuous strip chart of the gross gamma-count rate and 
the radar altitude as a function of time. The data collecting system 
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is shouTi in figure 3. An air filter sampler v;as also available on the 
aircraft to measure the concentration of radioactive airborne radon 
daughters which are the principal background contributors to the measured 
terrestrial radiation signals. 

Tedmique 

The detection system together with associated electronics was flown 
on a Beechcraft Twin Bonanza aircraft at an altitude of 150 m (500 ft) 
above the terrain and at a speed of about 140 knots. The two flight 
lines were surveyed three times in succession just after sunrise on 
March 18 and three times again in the early afternoon of the same day. 

The position of the aircraft and its altitude above terrain were measured 
and recorded every other second by a radar navigation computer system. 

The position and altitude measurements are accurate to +_106 m (j^350 ft) 
and m (+_20 ft), respectively. At an altitude of 150 m (500 ft) the 
field of view of the detectors covered an area of approximately 350 hectares 

for the mean energy of gamma radiation from naturally occurri.tg 
radioelements. 

All flights were timed to coincide with the NASA P-3 aircraft 
survey and the ground sampling work. The gross gamma counts were 
measured by the array of fourteen 4-in by 4-in Nal crystals nn board the 
survey aircraft. The data were recorded once every two seconds, pro- 
viding a data point for every 150 m (500 feet) of survey line traveled.' These 
data were averaged on a mile-by-mile basis. In addition, the gamma 
energy spectrum was recorded at eight-mile intervals during the survey, 
or twice per survey line. The spectral data were used to sura gamma 
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counts in the photopeak region and in the photopeak region. 

Soil moisture determinations were made using the three independent 
analyses of gamma flux attenuation from gross counts between SOkeV and 
3.0 MeV, as well as using the 40j^ photopeak gammas. All 

three sets of data, including previous measurements, were later employed 
in determining soil moisture values. 

In order to establish the magnitude of the airborne radon daughter 
gamma background during this survey, air filter samples uexe taken 
twice during the morning survey and twice during the afternoon survey. 

The filters trapped radioactive radon daughters which commonly attach 
themselves to aerosol particles in the size range between a fraction of 
a micron to a few tens of microns. Additionally, the gamma ray count 
rate was measured while flying the aircraft at 150 m (500 ft) over nearby 
Lake George where one would expect no terrestiial gamma radiation to be 
present--as another evaluation of the airborne radon daughter contribu- 
tion to the measurements. 

Results 

The gamma count rates in the photopeaks for and ^^®T1 were 
determined by first subtracting the non-terrestrial background counts 
(due to the aircraft, cosmic rays and airborne rd.don daughters) from 
the respective pulse height-spectrum windows (see figure 4) , i.e. , channels 
91 through 104 for and channels 168 through 184 for 208jj^^ These 
net photopeak count rates were then adjusted to correspond to that for 
transport of these terrestrial gajnmas through an air mass of 17 g/cm^ 
(equivalent to an altitude of 150 ra (500 ft) at average pressure and 
temperature). Table 2 gives a listing Of net photopeak count rates for 
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Relative Counts 


,40K Photopeak 


S 08 T 1 Photo peak 


Channel No. 


Energy (MeV) 


O I rHl Q ^ 

CO t- CO 

CO to CO 


C- ICO 
CO ICO CO 


to |> 


Definition of the gamma energy spectral pulse height windows 
used in the photopeak area equations. The letter labels ^ 
refer to the following: T (counts in the ^OS^i window); U. 

(counts in the i^indow occupied by the higher enei'gy 214Bi 
gammas from the uranium decay chain); B (counts in the 1.6 
Me V214si photopeak window) ; K (counts in the photopeak 
windoif) ; and L (lower energy background gamma counts around 
the 1.12 MeV^^'^Bi photopeak, (from Feimster et al, 1975) 


and for each survey line flown on March 18, 1975, and values 

averaged over repeated flights. 

The ratio of gamma flux rj/r 2 measured for two different soil moisture 
conditions Mj and M 2 has the form 

£1 = 1 + M2 Eq. (1) 

12 1 vMi 

Eq. (1) is used with the tabulated and ^^%1 photopeak count rates 
given in Table 2, the and photopeak count rate for November 

21, 1972, (2740 cpm and 583 cpm, respectively) and the soil moisture for 
November 21, 1972, (19%) to determine the soil moisture values for the 
March 18, 1975, survey. It should be noted that this can only be done 
for survey line 1 (west of 91st Avenue) since only that line was surveyed 
previously, a necessary prerequisite for determining the value of Mp 
of Eq. (1) . 

40 9ns 

The calculated soil moisture values, employing the K and the ‘‘^“Tl 
photopeak area methods are: 

M = 22% 

M (208ti) =20% ■ 

The third method of determining soil moisture from aerial measure- 
ments of terrestrial gamma radiation employs analysis of the gross gamma 
count rates between 50 keV and 3.0 MeV. If there is negligible airborne 
radon daughter contribution to the signal, the net terrestrial gross 
count rate is determined by subtracting the gamma background count rate 
due to the aircraft and cosmic rays. The net count rate is then adjusted 
to correspond to that for transport of 17 g/cm^ as was done for the 40j(^ 
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and 208^1 count rates. Occasionally, however, substantial airborne radon 
daughters are present during the surveys. Two approaches were used to 
establish the magnitude of the radon daughter background. The first was 
to determine this background contribution by utilizing (1) the change in 
the rario U/K of the gamma count rates in the two pulse height windows 
U and K (see figure 5), on November 21, 1972, and March 18, 1975, (2) the 
Ccounts/sec) /(counts/inin) due to airborne radon daughters in the U pulse height 
window as determined by experimental lake overflights, (3) a knowledge 
of the photopeak and radon daughter counts on November 21, 1972, and 
(4) a knowledge of the photopeak count on March 18, 1975. These 
values are, respectively 

1) U/K CNovember 21, 1972} = .44 U/K (March IS, 1975)'= 0.46 
2} Airborne radon count over water in U pulse height = 0.82 (counts/ 

sec}/(counts/rain) 

3} ^®K photopeak count (November 21, 1972} = 2740 cpm 
radon daughter count (November 21, 1972} = 600 cps 
4} ‘^^K photopeak count (March 18, 1975) = 2665 cpm 
With this data it is possible to determine the total average nonterrestrial 
gross count rate by proceeding through the following computations: 

Uranium terrestrial gross count rate (Tjovember 21, 1972} 

= 0.82 (U/K) (“^^K} - airborne radon daughter count 
count 

= 0.82 (0.44) (2740} - 600 (Eq. 2) 

= 389 cps 

)F THE 
POOR 
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TOLLESON 


TOLLESON 


91 AVENUE 


Portion of Landsat-1 image (3/17/75, 1967-17154, orbit 
13,484) covering Phoenix and the area to the west of the 
city including the test site q'^Tolleson, Arizona. MSS band 
5 (0.6-0.7pra) on left, MSS band 7 (0.8-1. 1pm) on right. 
Computer- generated and enhanced. 











Total uranium daughter (radon + terrestrial^ count rate (March 18, 1975) 
= 0.82 (U/K)(^®k) 

=0.82 (0.46) (2665) (Eq. 3) 

= 1005 cps 

gamma count rate (March IS, 1975) 

= (Eq. 2) 

- (Percent change in photopeak count between dates) 

X (Eq. 2) 


PFPRODUCIBILITY OF THE 
= 389 (0.03 X 389) PAGE IS POOR 


(Eq. 4) 


= 377 cps 

Radon gamma count rate (March 18, 1975) 

= (Eq. 3) - (Eq. 4) 

= 1005-377 (Eq. 5) 

= 628 cps 

Therefore, total average nonterrestrial gross count component for flight 
line is (radon + aircraft background + cosmic) 

=((Eq. 5) +290 + (137e(15-pxessure)jy2.55 (Eq. 6) 

= 1129 cps (at 13;90 lbs/in2barometric pressure) 

By determining the total 214gi photopeak areas in the U pulse height 
window for each flight line segment, approximate values for the radon and 
the total nonterrestrial backgroimd component for each line were computed 
as s]i own above. The results are tabulated in Table 3. The second method 
used to establish the contribution of the gamma count rate from the 
airborne radon was to compute an adjusted gross count (corrected for the 
nonterrestrial component) that corresponds to the neiv soil moisture 
computed for March 18, 1975, as described earlier. The adjusted gross 
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TABLE 3. Nonterrestrial background determined from the estimated radon 
activity using the U/K method. (Non terrestrial bkg = radon + 
aircraft + cosmic**). (From Feimster, et al., 1976) 


Gamma Counts* 
in the U pulse 
height windows 
(cpra) 

Pressure** 

(lb/in2) 

Radon Activity 
in the gross 
count channel 
(cps) 

Nonterrestrial bkg 
in the gross 
count channel 
(cps). 

938 

13.16 

556 

1151 

1027 

14.12 

645 

1222 

1223 

14.12 

S41 

1329 

1031 

14.12 

649 

1138 

1334 

14.16 

952 

1437 

1203 

14.16 

821 

1306 

1033 

14.09 

651 

1139 

1560 

14.13 

1178 

1666 

1376 

14.16 

994 

1522 

986 

15.64 

604 

1102 

* These data represent Lake ^ 

Roosevelt; Lines l.AN, 

2AS, IBN, IBS, 2BS, 

ICN, ICS, 2CS 

j and finally 

Lake Roosevelt. The 

line indication such 

as IAN means 

line 1, first 

pass, north half, etc 



**Cosmic radiation flux C is dependent upon the pressure according to 
C = 137 e (15-P)/2.55^ 


TABLE 4. Nonterrestrial backgi'ound computed from the adjusted gross 
count rate, fitted to photopeak area data. (From Feimster, 
et al., 1976) 

Fitted Adjusted Resultant 


Line 

Gross Count 
Rate (cps) 

Gross Count 
Rate (CPS) 

Nonterrestrial bk 
(cps) 

lA (AM) 

3682 

2500 

1307 

IB (AM) 

3780 

2500 

1399 

1C (AM) 

3776 

2500 

1389 

lA (PM) 

3837 

2500 

1381 

IB (PM) 

3526 

2500 

1070 

1C (PM) 

3065 

2500 

1009 




count rate for March 18, 1975, can be solved for by the equation: 


T2 = rjCl + Ml) 
1 + M 2 

where 


(Eq. 7) 


Tj and Mi = adjusted gross count rate and soil moisture, respectively, 
on November 21, 1972 

P 2 and M 2 = Unknown adjusted gross count rate and knovvTi soil moisture 
from photopeak area data on March 18, 1975. 

The average adjusted gross count rate and soil moisture during November 
21, 1972 was 2553 cps and 19%, respectively. The average soil moisture 
value determined from the photopeak area analysis of the March 18, 1975, 
survey data was computed to be 21.5% (average of tlie moisture values 
computed from the and 20Si>i photopeak areas) . Substituting tlie 
appropriate values into (Eq. 7), the adjusted gross count rate is: 


T2 = 2533 cps (1.190) = 2500 cps 
(1.215) 


(Eq. 8) 


To solve for the nonterrestrial component (radon + aircraft background + 
cosmic), the adjusted gross count is subtracted from the total gross 
count. The subtraction is accomplished by making use of the following 
two equations: 


GCjyjgT “ ^^TOT ~ (aircraft background + radon + cosmic) 
GC^ds = GC^et eG-054 (AM-17) 

where 


(Eq. 9) 
(Eq. 10) 


GGnet “ total gross count less the nonterrestrial component 
GCjoT ” total gross count 

GCajD “ gross count adjusted for air mass attenuation 
AM = air mass 

B-20 


,) is measured, and the expected adjusted gross 


The total gross count 
count (GC^qj) was established from the peak area soil moisture data 
(E.;i. . 8) . Combining equations gives : 

^*^ADJ “ (GC'pot ~ aircraft backgroiaid + radon + cosmic/2.55) 

X e0-054CAM-17) (Eq. 11) 

Solving for the nonterrestrial component leads to: 

Caixci'aft background + radon + cosmic) 

= - GC^pj e - 0-054 (AM-17) 12) 

Substituting the values for tlie aircraft and cosmic, the equation becomes: 


(290-+ Radon + 137o = GC^ - GC^pj e“0-54CAM-17) (Eq. 13) 


This allows computation of the total nonterrestrial component as well 
as the radon singly. Table 4 lists the values that are computed from 
this equation. Example 1 considers the first case in Table 4. 

Ex. (1) 

GCfOT ~ 4682 
AM =17.91 
GC^ = 2500 

(aircraft background + radon + cosmic) = GCiiQi-GC^pje" 0.054(AM-17) 

= 3682-2500e" 0.054(1791-17) 

= 1302 (nonterrestrial 
component) 


Substituting: aircraft background— 290 

cosmic = 137 e C15-p)/2.55 


Radon = 1302-290-137 e (15-14. 12)/2. 55 
= 819 cps 

This adjusted gross count rate was used for the tv/o-segment average of 
each pass made over the survey line listed in Table 4. 

Soil Moisture Ground Measurements 

Four points were sampled to a depth of 15 cm in each of 46 fields 
(16 vegetated, 30 bare) on each side of two flight lines in the following 
intervals; 0-1 cm, 1-2 cm, 2-5 cm, 5-9 cm and 9-15 cm. At each point, 
five samples were taken from the top and five samples from the bottom of 
the furrows. Average moisture contents for each of tJie four points were 
computed from each respective five interval samples. Soil moisture 
computed from soil' samples and weighted to a depth of 9 cm were chosen 
to compare with the aerial data. The weighing equation for the 9-cm depth 
is written as: 

% moisture (0-9 cm) = 1/9 (% mois. (0-1 cm) + % mois. (1-2 cm) 

+ 3 X % mois. (2-5 cm) + 4 x % mois. (5-9 cm) 

Table 5 lists all the areas that samples were taken from and their 
corresponding moisture values weighted for a 9-cm depth using the above 
relationship. 

The results in Table 6 show that moisture values computed from the 
gross count rate of the aerial data of March 18, 1975, are slightly 
higher than the ground-based measurements during the same period as 
measurements taken on November 21, 1972, and March 18, 1975. The average 
of the ground sampling moisture values was very nearly the same as the 
moisture value computed from the photopeak area data. 

As explained earlier, the gross count data are affected by the 


TABLE 5. Soil moisture values from soil samples weighted for the 
sampled 9-cm soil depth. (From Feimster et al., 1976) 


Average soil moisture values (% by weight of soil) 


Field 


North Point 


South Point 

97 (A§B) 


23.19 



22.53 

99 " 


19.78 



19.43 

101 »' 


19.78 



17.88 

102 " 


29.98 



31.75 

103 " 


23.59 



23.97 

104 " 


27.25 



28.98 

105 ” 


25.89 



26.81 

106 " 


30.30 



27.90 

107 " 


31.30 



53.78 

111 " 


22.89 



28.27 

112 " 


27.71 



29.74 

113 " 


27.27 



28.98 

114 " 


20,53 



18. ,69 

115 " 


24 . 77 



27.75 

116 " 


22.95 



21.22 

121 " 


15.56 



20.99 

123 f, 124 


24 . 10 



23.09 

126 (AfiB) 


19.17 



18.77 

127 § 128 


19.05 



19.40 

129 § 130 


16.90 



13.49 

131 § 132 


21.09 



19.94 

133 (A5B) 


9.38 



8.20 

134 ?T 135 


18.11 



14.47 

TABLE 6. Comparison of soil moisture values 

determined from ground 

sampling and 

from aerial 

measurements . 

(From Feimster, et al.. 



Soil 

Moisture (% 

by 

weight) 

Location 

Results 

from 

Results from 

Results from 

on line 1 

soil samples 

soil samples 

aerial measurements 

(mile no.) 

Nov 21, 

1972 

Mar 18, 1975 

Mar IS , 1975 

9 

21.3 


21 . 3 


16.6 

11 

28.8 


26.7 


32.2 

12 

22.3 


22.1 


25.6 

13 

18.2 


16,2 


19.1 

Avg. value for 

19.0 


21.56M.3 


24.3+6.1* 

entire survey 





21.5** 

line 






* These results 

are from the mile-by -mile gross 

gamma count rate analysis 


**These results are from the analysii; of gamma photopeak data, averaged 
over the entire survey line. 


airborne radon activity which is the major cause of the differences 
sho'.vn in Table 6 of the ground-based moisture determination and those 
from the aerial gross count data. 

Lands at Investigation - Visible and Near- I R 
Introduction 

Findings of previous investigators (Sabatini, et al, 1971; Steiner 
and Guterraan, 1966; Merritt and Hall, 1973; Planet, 1970; and Angstrom, 
1925) indicate that the reflectance of soils decreases in the visible 
and near-IR as soil moisture increases above the wilting point. The 
change in reflectance is also dependent to a significant degree on the 
soil type, with dark soils showing a greater change than light soils. 

The availability of Landsat-1 data (Orbit 13,484 ID 1967-17154) on the 
day preceding the ground observations at Phoenix permitted a test of 
the above-mentioned investigators' findings. 

Procedure 

The March 17, 1975, Lands at data v;ere displayed pictorially in 

Digital Muirhead Display images (DMD's) from 7-track CCT's with bands 

5 and 7 side-by-side (see figure 5). The test site along and on either 

side of 91st Avenue was located in the imagery and enlarged, figure 6a 

and 6b, to better view the area. A microfilm composite of the area in 

rectified 32 by 32 pixel arrays provided the quantitative data used in i 

the regression analyses. One such array is presented in figure 7 shoviring 
94, 95, 96a, 97A, 97B, 98A, 99A, 99B, and lOOA. 
fields / Even though the images appear "blocky" due to the enlargement, 

the 80-meter resolution of the MSS still allowed 25-35 satellite data 

points per field. Variation in the number of points was a function of 

a given field's boundaries with respect to the scan lines of the MSS. 
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TOLIESON 


AVENUE 








Enlarged segment o£ MSS band 7 shown in Figure 5, covering 
test site and the immediate vicinity around Tolleson, Arizona. 


Fig. 6b. 
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The satellite values were averaged for each field for direct 
comparison with the averaged soil moisture values taken from Blanchard's 
(1975) report. As the reflectance measured by the '-BS is essentially 
from the surface of the vegetation or soil, field soil moisture values 
of the top one centimeter (for both top and bottom of the furrow) from 
the following morning's data survey (March IS, 1975) v;ere used to corre- 
late with Landsat MSS readings. The difference in time between the two 
data sets--about 24 hours — is not as critical as comparing early morning 
and late afternoon values on the same day. Jackson et al. (1975) showed 
that there is a recovery of soil moisture , an upward flux of moisture, 
beginning several ]iours before sunset, to partially replace those losses 
(up to 75% v/hen the soil surface is relatively dry) due to evapotranspifa 
tion during the day. Tlie recovery may be seen in figure 8. Satellite 
MSS reflectance values and soil moisture values for vegetated and bare 
fields were compared . 

Results and Discussion 

Plots of the various MSS band reflectances versus soil moisture in 
percent are presented in figures 9a- d and 10. There appears to be little 
change in reflectance for vegetated fields in either of the visible bands 
MSS 4 and 5. However, in the near-IR (MSS 6 and 7) , there is an increas- 
ing reflectance in both bands with increasing soil moisture for alfalfa, 
but not for I wheat . Apparently the additional moisture is manifest in a 
more vigorously growing alfalfa plant, giving a brighter response in- 
bands 6 and 7. The same does not seem to be true for wheat. 

All bare fields showed a decrease in reflectance as soil moisture 
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Fig. 8. Volumetric water content in the 0.0 to 0.5 -cm increment 

versus time for 3 days during March 1971. The solid line 
represents smoothed data and the symbols represent the 
measured values, (from Jackson, 1973) 
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Fig. 9b. Regression analysis of percent soil moisture v,s. MSS reflect- 


ance raw counts, MSS band 4, vegetated fields only 
(parabolic fit). 
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PERCENT SOIL MOISTURE, BY WEIGHT 


Fig. 9d. Regression analysis of percent soil moisture vs. MSS reflect- 
ance raw counts, MSS band 6, vegetated fields only, alfalfa 
(parabolic fit) and wheat (power fit) separated. 








RELATIVE REFLECTANCE. ARBITRARY UNITS 



Fig. 10. Combined plots of MSS bands 5 and 7 raw counts vs. soil 

moisture: a) MSS band 5, bare fields only (power fit); 

b) MSS band 5, vegetated fields only (straight-line fit);- 

c) MSS band 7, bare fields only (power fit); d) MSS band 7, 

j ^ alfalfa only (par'’’^olie fit); and e)MSS band 7, wheat only 

If fit). 



increased, with a power fit to the data consistently supplying the best 
correlation of the five available curves for data fitting. The amount 
of variance, r^, ranged from 0.71 for MSS 4 to 0.76 for MSS 6 and 7; r^ for 
MSS 5 was 0.75. Even though the fits were good, most of the reflectance 
decrease occurred between 5 and 20 percent soil moisture with very 
little further decrease from 20 to 55 percent. 

Only 18 data points were available for vegetated fields, 8 from 
fields growing alfalfa and 10 from fields in wheat. The two visible 
bands in the MSS, 4 and 5, showed n£ change in reflectance as soil mois- 
ture varied from 5 to 55 percent. The brightness of wheat fields in the near-IR 
(MSS 6 and 7) also indicated almost no dependency on soil moisture changes. 

The alfalfa fields in the near-IR increased in reflectance as the soil 
moisture increased. Unlike the curve for bare fields, near-IR reflectance 
values did not level off at high soil moisture. Rather the reflectance 
increased at a slightly increasing rate. One possible explanation for 
the differences noted between wheat and alfalfa is that the alfalfa 
provides a solid canopy; but the wheat, being only 22-50 cm (9-12 inches) 
high, did not provide a complete vegetative cover, and thus the satel- 
lite was integrating the response of both bare soil and vegetation. The 
alfalfa apparently responded to the increased soil moisture by growing 
more vigorously and appearing more reflective. 

Though visible and near-IR remote sensors provide some indication of 
soil moisture at the surface of bare fields, the correlations are best 
between 5 and 20 percent soil moisture. 

Figure 11 shows Landsat MSS band 7 raw data (full resolution) 
superimposed on a false color near-IR aerial photograph. The aircraft 



Fig. 11. False-color aerial photogaaph of fields 94, 95, 97A, 97B 
as well as portions of fields 91, 92, 93A, 98A, 99A, 99B, 
and lOOA, overprinted with MSS band 7 raw data from 
Landsat-1 image (3/17/75, 1967-17154). See Fig. 7. 
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data ivere collected on the afternoon of March 18 near 1 P.M. local 

time; the Landsat values were obtained at 10:00 A.M. local time on 

March 17. The area is west of Phoenix and covers a portion of the 

test site located north of Tolleson. 91st Avenue runs vertically (N-S) 

through the right (east) portion of the image. Red areas indicate 

actively growing vegetation. The large rectangular area just below the 

center of the image is test fields 97A and planted in alfalfa. Acquisition 

of soil moisture values, obtained close to the time of aircraft overpass, differed 

from the Landsat-acquired near-TR information by approximately 27 hours. 

However, as shown by JacKson et al. (1975), the effects of differences 
in observation times are somewhat compensated by the nocturnal 
replenishment of surface soil moisture from reserves deeper in the 
ground. Surface (0-2 cm) moisture values (in percent of dry weight) for 
fields 97A and B ranged from 13.70 to 27.65. The driest areas occurred 
in field 97A adjacent to 91st Avenue; the wettest along .the western 
(left) edge of field 97B. As shown in figure 10 a general increase in 
Landsat MSS band 7 values for alfalfa fields occurs as soil moisture increases. 

MSS band 7 values range from 64 in the drier portion of field 97A to 84 at i 
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jthe wettest locations in .field' 97B. This proportional relatipnshi£_ of soil 
moisture to MSS band 7 values was found to exist for all alfalfa fields 
sampled in the study (see figure 10) . 

The amount of useful information from visible and near-IR is 
miniscule when compared to that provided by passive microwaves and 
potentially from the thermal IR. Further efforts to predict soil 
from remote platforms should place the emphasis on the microwave portion 
of the spectrum and thermal IR where the most favorable results have 
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been achieved. For vegetated fields, though the MSS data were limited, 
crops growing complete canopies such as does alfalfa, may manifest 
increases in soil moisture by growing more vigorously. 

NOAA Investigation - Thermal IR 

Another area of the spectrum that might provide some information on 
soil moisture and its change, w’ith time is the thermal IR at 10.5-12.5vira. 
Work done by Idso, et al. (1975) indicate that if the soil type is known, 
the amplitude of the diurnal surface soil temperature wave is a function 
of the water content in the surface layer. Maximum values of surface 
soil tempei'atui'es were observed about 1.5 hours after local noon; minimum 
surface soil temperatures ju.st before sunrise. 

The 30-minute interval betv;een observations by MOAA's SHS/GOES satel- 
lites allows measurements of temperatures near the daily extremes and 
thus could closely measure the daily amplitude of surface soil temperatures 
A computer program, has been developed by J. Pritchard (personal communica- 
tion) that subtracts any two SMS /GOES tapes. An example is shown in 
figure 12. The difference between the two indicated observations is in 
the center, the subtraliend on the left, and the minuend on the right. 

Unfortunately, two problems exist with the data that preclude its 
usefulness as a potential soil-moisture monitor. Figure 13 vividly 
displays the problem of low resolution. In this image the blockiness 
of the individual 8-km resolution pixels is evident. Moreover the 
entire Phoenix test site area does not even fill one 8 by 8 km pixel 
(being only 3 km by 8 km). The test site's approximate location is 
indicated, as are other more easily identified landmarks. The second 
problem involves the difficulty in obtaining taped SMS/GOES data. This 
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SMS-2 8-km resolution therraal-IR (10.5-12.5 m) iciagery of 
southwestern North America and eastern Pacific Ocean. Image 
on left was observed at 2015 GMT on 3/17/75. Image on right 
was observed at 1015 GMT on 3/18/75. Image in center was 
obtained by subtracting temperatures measured on 3/18/75 
from those measured on 3/17/75. 
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Fig. 13. Enlarged portion of 8-km therraal-IR (10.5-12.5 in) image 
acquired by SMS-2 at 1015 GNfT on 3/18/75. Area viewed 

covers southwestern U.S. (including Phoenix test site) and | 

portions of Baja, California^ and the Gulf of Baja. ^ 
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example was obtained only by special request, as GOES tapes are not 
archived. The reason for not saving tapes is lack of storage facilities. 
Thus until tapes at an improved resolution from thermal sensors on 
geosynchronous satellites become a reality, the potential for a remote 
thermal soil moisture monitor will not be realized. 

Conclusions 

Estimates of soil moisture were obtained from visible, near-IR 
gamma ray and microwave data. Attempts using GOES the rmal-IR were 
unsuccessful due to poor resolution (8 km) , even though laboratory 
studies support the thermal inertia concept as a valid means to measure 
soil moisture fluctuations, Microwaves were the most effective at soil 
moisture estimates, with and without vegetative cover. Gamma rays 
provided only one value for the entire test site, produced by many data 
points obtained from overlapping 150-meter diameter circles. Even though 
the resulting averaged value was near the averaged field moisture value, 
this method suffers from atmospheric contaminants (randon daughters), the 
need to fly at low altitudes (150 ra) and the necessity of prior 
calibration of a given site. Visible and near-IR relationships are 
present for bare fields but appear to be limited to soil moisture levels 
between 5 and 20 percent. The densely vegetated alfalfa fields 
correlated only with near-IR reflectance. Soil moisture values from 
wheat fields showed no relation to either visible or near-IR MSS data. 


Summary of soil nioisttire 


The study, involvir4g aircraft and remote s,ensors, evaluated the 
effectiveness of. portions of the electromagnetic spectrum to estimate 
soil moisture. Tests were made at tw sites, viz., Phoenix Arizona, 
and Luveme, Minnesota, using data collected by visible, near-IR, 
thermal IR, microwaves and gamma-ray sensors. The results of these 
studies indicate the greatest potential for remote soil moisture monitors 
to be thermal IR and microwave sensors. Operationally both suffer from 
a lack of resolution. Microwaves can also penetrate vegetation and the 
soil surface to provide average soil moisture estimates to almost 8 
centimeters deep. Thermal inertia estimates are of necessity obtained 
from bare fields and can only indireGtly determine moisture conditions 
of the subsurface layer. 

Findings of the gamma-ray surveys are within range of field data. 
However, the method has no future in any satellite program, b^eihg limited 
to very low altitude (150 m) aircraft flights. The technique necessi- 
tates at least one previous flight over a given region to calibrate the 
system. Noise from radon daughters must be taken into account. The 
gamma-ray soil moisture technique does provide a check for systems such 
as remote passiye microwaves where coarser resolution prevails. 

Visible MSS data from Landsat correlated with soil moisture | 

only when the soil was bare, and vvould thus be a very limited tool to 
maintain year-round monitoring of soil-moisture variations. The near-IR 
appears someitfhat more promising showing correlation with bare ground and 
some Correlation with changes in soil moisture when a vegetative cover 
is present (alfalfa) .But results for wheat, for sparse data 


available, indicated no relationship between satellite near-IR reflect- 
ance and soil moisture. The overall results of the Landsat MSS data are 
marginal at best and are far inferior to those soil-moisture estimates 
obtained by the microwave technique. 

Future 

The continued development of sensors for future satellites iprovides / 
the hope tha^ a better understanding of soil moisture variations will 
be possible via remote sensing. In particular there were five satellites 
launched in 1978, each containing one or more sensors collecting data 
pertinent for soil moisture evaluation. Landsat-5, launched on March 
5, unlike its predecessors, has a thermal channel in the 10.5-12.5ijm 
portion of the spectrum. Unfortunately outgassing problems have 
continually plagued the thermal channel resulting in an increased NEAT 
(noise equivalent temperature increment) now- exceeding 2.50C. This 
degradation will, at some point, make the data totally useless, likely 
by June 1979. 

Though still 1-km in resolution, the daily coverage provided by the 
thermal-*IR of the AVHRR (Advanced VHRR) on TIROS-N, placed in orbit on 
13 October 1978 has a lower noise levei, only 0.1°C NEAT (noise equivalent 
difference in temperature) , than the IR scanners on board earlier NOAA , 
satellites. The HCMR (Heat Capacity Mapping Radioineter) on board the 
HCMJ'l satellite, oribiting the Earth since 26 April 1978, has a 500 ra 
resolution with a Ih to 3%-day cycle (8-day cycle for thermal inertia measure 
ments) . The thermal information collected by these three satellites 
will permit the furtherance of studies such as the one presently being 
done in Phoenix (Jackson, et al., 1975) , The HCMM satellite scattered 


a significant loss of power on 28 February 1979 resulting in restricted 
data collection beginning 10 March 1979. 

The other two satelliies^ launched in 1978 are Seasat-A and Nimbus- 
G. Both satellites collected inforraation of Earth in the microwave 
region of the spectrum. Seasat-A was launched on 26 June 1979; Nimbus-G 
on 24 October 1978. 

Both Nirabus-6 and Seasat-A have onboard the Scanning Multi- frequency 
Microwave Radiometer (SMMR) . Unfortunately the resolution (25 to 125 
km) and wavelengths (6.6, 10. 69, 18, 22.235, 36 GHz) are not optimum 
for soil moistrre evaluation. The active microwave system on Seasat-A, 
known as SAR ._.CS)aiethic ..Aperture Radar) imaged ocean and limited land 
areas at 1.35 GHz with a resolution of 25 m. This very high resolution 
imagery, though not in the optimum soil moisture wavelength range of 
4.25 to 5.25 GHz, should provide useful information regarding soil 
moisture evaluation, particularly at areas the size of agricultural 
units (65 hectares) . A massive power failure occurred 10 October 1978 
terminating all functions of the spacecraft. 
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.LUVERNE, MINNESOTA SOIL MOISTURE SURVEY 
Description of Test Site 

The Luverne soil moisture survey area lies in the extreme southwest 
corner of Minnesota, approxiTOately 30 miles northeast of Sioux Falls, 
South Dakota (see figure 14) , The test site is 8 miles long and is south 
of Luverne, Minnesota, and parallel to U.S. Highway 75 (see figure 15). 
The elevation of the test area is between 1,400 and 1,500 feet above 
mean sea level with slopes varying from 0 to 8%. Soil types consist 
mainly of silty clay loams and alluvial soils located near flood plains 
(see figure 16) and most of the acreage grown in corn, alfalfa, hay, and 
wheat (USDA, 1949). 

:i>oil Moisture Experiments - 

During the pei'iod of this study two soil moisture experiments were 
attempted in the Luverne test site area. The firstwas performed May 
11-12, 1975, and the second on October 17, 1975. During both of these 
test periods, ground truth, meteorological conditions and aircraft and 
satellite remote sensing platforms were to be used to evaluate soil 
moisture conditions. Unfortunately during the May 11-12 test period 
cloud contamination prevented the use of satellite data and during the 
October 17 study no Landsat overpass was scheduled and no ground truth 
was collected. 

The May 11-12, 1975 Experiment 

On May IT, light rain showers persisted until approximately 9:30 a, iji. 
After these early morning showers, slow clearing commenced and by noon 
mostly clear skies prevailed. By 2 : 00 p .ra. , however) scattered showers ' 
and thunderstorms developed and continued throughout the afternoon. 
National Weather Service records for Luverne indicate 0.47 inches of rain 
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Fig. 16. Soil map of Luveme, Minnesotaitest site. U.S. Highway #75 
is shown by dark bold line. Test site line is shown by 
dashes. Map taken from USDA soil map of Rock County, 
Minnesota. 
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for the day. Because of the inclement weather the planned EG§G gamma- 
ray aircraft underflight was cancelled, the Landsat overpass imagery 
proved unusable dub to cloud cover over the test site, and no ground 
truth was collected. 

Fortunately May 12 was without inclement weather and the EG$G aerial 
survey for the measurement of natural terrestrial gamma radiation was 
carried out. The purpose of the survey was similar to the passive 
gamma ray survey performed over the Phoenix test site, i.e., to determine 
soil moisture by detecting differences in the terrestrial gamma flux 
attenuation as a function of soil moisture content. Initially it was 
hoped to establish correlations betvveen 1) the NASA Landsat-2 satellite 
measurements of the terrestrial near-infrared and optical signals and 
2) the soil moisture determined from ground-based sampling. The EG§G 
aerial measurements of terrestrial gamma radiation were planned to be 
simultaneous with the NASA satellite measurements and the ground 
sampling thus providing an independent set of data from which to establish 
correlations and evaluate the NASA near-infrared data. The planned com- 
parisons were not possible, however, because, as mentioned previously, 
of heavy cloud cover during the time that the satellite was due to pass 
over the area. 

A. Equipment 

I The equipment used for the Luverne gamma ray survey was the same 

as that used for the Phoenix gamma ray survey, (see p. B-lO of this report 
for details) . 

B. Technique 

This section on technique is taken largely from Feimster and 
Fritzsche, 1975. The Luverne , Minnesota, survey included five lines 


approximately 8 miles long, three flovm in the north-south direction and 
two in the east-west direction, oriented parallel and perpendicular to 
U,S. Highway #75, respectively. Radiological data used to compute soil 
moistures were collected over Line A, shown in figures 17 and 18. This 
line is located 1000 ft west and parallel to Highway #75. 

A Twin Beech Bonanza aircraft carrying radiation detection instrumen- 
tation as sho\m in figure 5, and flying at an altitude of 500 ft (152 m) 
above the terrain, made four passes over each of the five survey lines. 

The flights were timed to coincide with the ground sampling work. 

Details of the technique are given in the Phoenix test site section 
and are not repeated here. 

C. Results 

The predominant gamma rays from natural radioelements in the 
soil include the 2. 6MeV 205-j-j^ gammas from the thorium decay chain, a 
family of ^^^Bi gammas from the uranium decay chain ranging over values 
of 0.61, 0.76, 0.93, 1.12, 1.76 and 2.2 MeV, and the I. 46 MeV gammas 
from '^Ok. a typical terrestrial gamma-ray pulse height spectrum as 
measured by a large Nal(Tl) crystal is illustrated in figure 4. The 
gamma counts under the photopeak, between pulse height channels numbef 

168 and 184> contain no background contribution from the airborne ^^%i | 

radon daughters, since the 214gi ganimas are at lower energy values; the 

photopeak data provide an abundance of counts, and can be separated 


from higher energy 214 bi corapton tails through a spectral stripping 


analysis. The gross count data, integrated over the range between 50 keV 
and 3.0 MeV offers the advantage of excellent counting statistics and 




Luveme, Minnesota^ survey flight pattexn-north half . The 
dashed lines labeled A through C and E indicate the positions of 
survey lines (from Feiraster and Fritzsche, 1975). 
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TABLE 7. Not 40 k and 20Sjj gamma photopeak count rates from the May 12, 1975 
survey of Line A at Luverne, Minnesota (from Feimster and 
Fritzsche, 1975) 


Spectrum 

No. 

Pass 

No. 

Net Gamma Count 
40 k 

Rate in Photopeak (cpra) 
20§xi 

55 

1 

1859 

556 

60 

2 

1907 

540 

65 

■ 

o 

. 1841 

537 

71 

4 

1905 

542 



Avg. 18/8 

Avg . 544 



photopeak area methods are: 

M C^^K) = 32% 

M C^O^Tl) = 28% 

The third method of determining soil moisture from aerial measurements 
of terrestrial gamaa radiation employs and analysis of the gross gamma- 
count rates betiveen 50 KeV and 3.0 MeV as explained in the Phoenix 
test site section. Correcting the gross gamma count for the airborne 
radon contribution is most often achieved by air filter data analysis 
(See Feimster and Fritzsche, 1975 Appendix A) . Table 8 lists data 
collected by filter and radon values computed from these filter data of 
May 12, 1975 and March 6, 1973. 

To validate the air filter data, a method called U/K ratio was 
employed. The significant elements of this method are: 

1. The uranium photopeak to potassium photopeak count rate ratio 
is computed for any two surveys under consideration. These are labeled 
U and K as parts of the gamma energy spectrum defined in figure 6. 

2. If the U/K ratio increased from one soil moisture survey to the 
next, then the radon daughter contribution to the gamma count rate has 
increased in' both the energy spectra and in the gross count rate. The 
U/K ratio is unaffected by soil moisture changes. 

Table 9 lists the U/K ratio computations for the March 6, 1973 and 
May 12, 1975 data. Since the U/K ratio decreased from the first to the 
second survey, then the airborne radon must have decreased. This resull 
supports the air filter results (in Table 8) which show the radon contri- 
bution decreasing from 600 cps to 165 cps. 

The air filter method is used fo^ the soil moisture analysis in this 
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TABLE 8. Air Filter Data (from Feimster and Fritzsche, 1975) 


Excess^ 


Date 

Air Filter 
Comt Rate 
(per 3 min.) 

Background 
Count Rate 
(per 3 rain.) 

Net Filter 
Count Rate 
(per 3 rain.) 

Gross Gamma^ 
Count Rate 
(sec“^) 

Filter 2 
Factor 

Location 

Gross Gamma 
Count Rate 
(sec~^) 

5/12/75 

4260 

962 

3298 

1100 

.18 

Lake 

165 


1794 

877 

916 



Land 


3/6/73 

3778 

1617 

2161 

1200 

.32 

Lake 



3524 

1604 

1920 



Land 

600 

^The gross 

gamma count 

rate is the gamma 

count rate per 

second in the 

primary detector array. 



^The filter factor is the gross gamma count rate less tlie aij'craft and cosniic ray contributions divided 
by the net filter count rate. 


^The excess gross gamma count rate is due to airborne radon daughters. It is computed from the filter 
factor times ; the net filter count rate. 


TABLE 9. Gamma spectral data used in U/K ratio computations (from Feimster 
and Fritzsche, 1975). 


Survey 

Date 

Gamma count rate in 
K (cpro) 

selected pulse height window* 
U Ccpm) 

Ratio 

U/K 

3/6/73 

1719 

1285 

0.75 

5/12/75 

1877 

616 

0.33 


♦Refer to pulse height window definition in Figure 4, 


report, because it yields an absolute magnitude for airborne radon. The 
U/K ratio yields the change in concentration for airborne radon only. 

The contributions of airborne radon to the primary detector signal as given 
in Table 2 are: 165 cps for May 12, 1975 and 600 cps for March 6, 1973. 

The cosmic + aircraft contribution was 500* cps. Therefore, the non- 
terrestrial background was 1100 cps on March 6, 1973 and 665 cps on 
May 12, 1975. These non-terrestrial background count rates were subtracted 
from the gross count rates along each mile of survey line, giving the 
net terrestrial gi'oss count rate averaged over a Survey distance of one 
mile. Table 10 shows the adjusted gross count rate per mile for four 
passes made over Line A, The adjusted gross count rates averaged raile-by- 
raile basis as shorn in Table 11. 

To support the soil moisture values computed from aerial data, 81 

soil samples were taken and their moisture content measured. Moisture 

values from these samples were supplied by Twin Cities Testing and 

Engineering, Sioux Falls, South Dakota. Their method was to weigh the 

samples before and after drying and compute the soil moisture values 

determined from the ratio of the wet minus dry weights divided by dry 

weights. Table 12 compares the moisture values on a mile-by-mile basis 
(Table 12a) 

and a 2-mile basis/i determined from the aerial data and the moisture 
values from the ground sampling data for the two survey dates. Some 
discrepancy is generally e:q)ected between the ground and aerial data. 

As explained earlier, moisture values acquired by the aerial system are 
computed form data averaged over approximately a % sq. ma. area opposed 
to point samples taken to compute ground based soil moistures. 

r As a check for data acquisition consistency of the aerial system, 

several passes were made over each survey line. A plot, figure 19, of 

*Cosmic count rate = 290 + 137e(15 - pressure) /2. 55) cps . 
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TABLE 10. Adjusted and averaged adjusted gross count rates per mile for 
line A (from Feimster and Fritzsche, 1975). 




May 12, 1975 


March 6, 1973 

Pass No. 
Mile 

1 

2 

3 

4 

Average 

Average 

1 

2293 

2315 

2359 

2262 

2307 

1930 

2 

2400 

2484 

2377 

2343 

2391 

1969 

3 

2553 

2499 

2531 

2501 

2508 

2195 

4 

2570 

2552 

2552 

2531 

2552 

2304 

5 

2492 

2530 

2486 

2485 

2498 

2248 

6 

2549 

2578 

2519 

2491 

2534 

2266 

7 

2443 

2496 

2447 

2413 

2450 

2308 

8 

2462 

2441 

2397 

2447 

2436 

2239 


Note: Units of gross count are counts sec"^. 


Table 11. Soil moisture for May 12, 1975 computed from averaged adjusted 
gross count rates in Table 10 and March 6, 1973 soil moisture 
data (from Feimster and Fritzsche, 1975). 


March 6, 1973 May 12, 1975 



Gross Count 

% Soil* 

Gross Coul^t 

% Soil^" 

Mile 

Rate (sec"l) 

Moisture 

Rate (sec"^) 

Moisture 

1 

1930 

47 

2307 

23 

2 

1969 

47 

2391 

21 

3 

2195 

46 

2508 

28 

4 

2304 

51 

2552 

36 

5 

2248 

41 

2498 

27 

6 

2266 

41 

2534 

26 

7 

2308 

43 

2450 

35 

8 

2239 

: 35^V ■; 

2436 

24 

Averages 2182 

44 

2460 

28 


*The soil moisture on March 6, 1973, was obtained from ground-based 
measurements. 


‘‘’The soil moisture on May 12, 1975, was obtained from the airborne 
gross count rate data from both March 6 and May 12 and the soil 
moisture data of March 6, 1973 (see Equation 14) . 

REPRODUCIBILITY OF THE 
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TABLE 12. CorapaTison of the aii'horne gross count rate method and the 

ground-based soil moisture measurements per mile (from Feiraster 
and Frittsche, 1975). 


Soil .Moisture (%) 

Mile .Airborne Ground Based 


1 

23 

25 

2 

21 

26 

3 

28 

39 

4 

36 

28 

5 

27 

23 

6 

26 

27 

7 

35 

30 

8 

24 

29 


Average 28 

28 

T.ABLE 12 a. 
Miles 

Comparison of airborne gross count rate 
ground-based soil moisture measurements 
Feimster and Fritzsche, 1975). 

Soil 

Airborne 

method and the 
per 2 miles (from 

Moisture 

Ground Based 

1-2 

22 

25.5 

3-4 

32 

35 

5-6 

26.5 

25 

7-8 

30.5 

29.5 


Average 28 

28 


Distance along survey line A - miles 


Average adjusted gross count rate per mile for each of the 
four paisses of line A (from Feimster and Fritzsche7~1975) . 




^ Q 

rO 
' 9 * ' 


-■feiti wjLifc 'Si ’ . V 


4 - 



the variation tlje adjusted gross gamma count rate during each pass over 
surve>^ Line A, using gross count rate versus time, illustrates acquisition 
consistency during this survey. The difference between corresponding 
points of each curve is due primarily to counting error which is 
approximately (count rate)^/ (count rate) at a 68% confidence level; 
i.e., + one standard deviation. 

Figure 20 shows the relationship between the ground sampled soil 
moisture values values computed per mile and the soil moistures per mile 
from the adjusted gross count aerial data. 

D. Conclusions 

The airborne gamma-ray results indicate a range of vSues for 
soil moisture for the sui'vey at Luverne, Minnesota, on May 12, 1975, from 
28% to 39%. The spread in these values is due to: 

1. The statistical counting uncertainty of the photopeak area 

methods. 

2. The airborne radon concentration uncertainty in the gross 
count rate method . 

, 3. The uncertainty in the ground me as ureraerits of the baseline 
soil moisture, i.e., that soil moisture obtained from soil samples on 
March 6, 1973. 

The and 208^1 pi^otopeak methods gave average soil moisture values 
of 32% and 28%, • respectively . The adjusted gross-count rate method .eave ! 
an average soil moisture value of 28%, with mile-by-mile values ranging 
from 21% to 36%. Ground measured soil moisture values ranged from 23% 
to 39% with the average value being 28%. Thus, the three methods of 
surveying soil moisture by airborne y-ray gives results consistent with 
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those obtained by ground truth. The photopeak methods, however, appears 
more useful in measuring the long-line soil moisture average, whereas the 
adjusted gross-count' rate method can be used to measure soil moisture in one 
mile increments. 

To improve the accuracy of the airborne gamma-ray soil-moisture measure- 
ments, the following experiments and procedures are recommended: 

a. Increase the amovint -j soil samples taken to reduce the true soil 
moisture uncertainty. 

b. Improve the systematic method of recovering soil samples with respect 
to location and sampling density on the survey line. This could be accom- 
plished by permanently marking the soil sample sites for future sampling. 
However, due to crop rotation and field tillage, this may be impractical except 
for selected site locations . 

c. Reduce the error in the airborne system due to radon. To accomplish 
this, the variable spatial distribution of radon daughters in the air would 
need to be determined, as \tfell as determining the varying relative humidity 
and temperature more accurately, i 

d. Establish a more accurate relationship of gamma-ray flux to soil 
moisture through more soil moisture surveys or ground-based experiments in a 
laboratory-controlled situation. 

The October 17, 1975 Experiment 

^ 1975, the NASA P-3 aircraft was scheduled to over- 

fly the Luverne, Minnesota, test site. The sensors aboard include an 
11-channel multispectral scanner, six AMPS cameras and filters, one Zeiss 
camera, and one PRT-5 radiation thermometer. The characteristics of the 
scanner and camera systems are listed in Tables 13, 14, and 15, respectively. 
The PRT-.5 radiation thermometer is used to measure gross apparent 
radiation in the 8 to 14pm range. The instrument has a 2° field of view 


TABLE 13. Characteristics 


Channel 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Spectral Bandwidth 
.370-. 413 
.440-. 490 
.495-.S35 
.540-. 580 
.580-. 620 
.620-. 660 
.660-. 700 
.700-. 740 
.760-. 860 
.970-1.060 
8.000-13.000 


Remarks 


Inoperative 


, ' 0 '" 

Inoperative 
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and a response time of 5 milliseconds. Incoming radiated energy, continu- 
ously compared to a 55°C internal reference, is converted into a voltage, 
which is directly related to the energy difference between the target 
and the reference. This voltage is displayed on the indicator as 


equivalent blackbody temperature to an accuracy of 0.5°C. 

Unfortunately, because of a very short notice that the P-3 aircraft 
was to overfly the test site area, there was no time to organize a ground 
data collection effort. Also no Landsat overpass was available for this 
day. Thus the resulting aircraft data cannot be verified by using ground 
data or satellite information as planned. 

Weather over the test site was clear and temperature at the time of 
the overflight was near 16°C. Weather records at nearby Sioux Falls, 

South Dakota, indicate little or no precipitation from October 14-17. 

Three overflights (lines 7, 8, and 9) were made over the test site 
parallel to U.S. Highway #75. Of these three passes line 7 had the test 
site line too near the edge of the imagery, line 8 has the test site 
line somewhat off center of the imagery, and line 9 was centered on the 
imagery. Owing; tp delays in securing data tapes and problems with the 
computer reduction of line 9, line 8 was analyzed instead of line 9. 

A. Zeiss and AMPS Data 

; ■ ; ; ' ■■ : ; I ^ 

Of the seven Zeiss and AMPS camera stations available as data 

sources on the October 17, 1975, flight, only four were ordered for 
analysis. The filter band passes selected were 400-700 nm (AMPS aerial 
color) , ^ 500-800 niti (AMPS color IR) , 875-925 nm (AMPS B§W IR) , and 
520-850 nm (Zeiss color IR), After examining the imagery, four fields 
labeled A, B,,C, and D, were selected for further study (figure 21). Thq 
fields, shown in the Zeiss and A^'IPS images (figures 22 through 39) were 
selected because 1) the imagery showed them to be devoid of vegetation, 
and 2) comparison of the soil types also showed them to be sirailar, 
i.e., loam. Specific soils for each field are listed in Table 16. 
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photograph of field A 
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AMPS aerial color photograph of field C 
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AMPS color IR photograph of field A 
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Fig. 31. AMPS color IR photograph of field B. 
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AMPS color IR photograph of southern section of field D 
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Fig. 36. MIPS black and white IR photograph of field B. 
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Fig. 37. AMPS black and white IR photograph of field C 
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Fig. 33. AMPS black and white IR photograph of the northern section of 
field D. 
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of the southern section of 


AMPS black and white IR photograph 
field D. 
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Since the four fields are similar, any difference between them can only 
be attributed to texture or soil moisture differences. The Zeiss and 
AMPS imagery, however, does not permit quantitative analysis of the 
field texture or soil-moisture conditions. Therefore, a quantitative 
analysis to detect texture or soil moisture differences, if any, between 
three fields was undertaken using the 11-channel multispectral scanner 
digital data. 

B. Multispectral Scanner Data 

Condit (1970) has shown through laboratory studies that general soil 
types and wet and dry soil conditions can be differentiated using five 
characteristic wavelengths; 0.440ym, 0.540pm, 0.640pra, 0.740pm, and 
0.860pm. From these five wavelengths a characteristic curve relating 
percent reflectance to wavelength can be constructed for three general 
soil types, i.e., chernozem, pedalfer, and laterite. Furthermore, 
these curves show certain changes when the soil type is changed from a 
dry to a wet state. Since the multispectral scanner onboard the NASA 
P-3 aircraft included the five wavelength channels used in Condit’s 
experiments, a study w'as performed to determine if the results could be 
duplicated using the four bare fields of the October 17 experiment over 
Luveme. According to Condit’s analysis the four bare fields should 
have a curve similar to that of figure 40, a chernozem type soil curve 
of which the soil types of the four fields are included. 

To derive the curves from the scanner data, the digital values for 
all wavelengths were converted to radiance values using an algorithm 
supplied by NASA/JSC. Contoured microfilm of the radiance values of 
the four fields was then generated for seven wavelengths of which the 
five used by Condit were included . The raid-wavelength values used were 
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0.465yni, 0.515vim, 0.560pm, 0.640pm, 0.680pm, 0.810pm, and 1.015pm. Two 
other visible spectral channels (0.600 and 0.720 m) were available on the 
M^S but were not used due to malfunctions in both channels, A typical con 
toured radiance value microfilm product is shown in figure 41. 

Plots of the average field radiance vs. wavelength for fields A, 

B, C, and D are shown in figures 42-45. Table 17 is a tabulation of the 
average field radiances for each field vs . wavelength as well as the 
mean and standard deviation of the radiances for a specified wavelength. 
Using Table 17, a comparison of the mean radiance values and standard 
deviations of each spectral bandwidth to the average I'adiance value of 
each field shows that only field D differs statistically from the 
other fields. Since the soil t>'pes of the four fields are the same, the 
possibilities accounting for this difference include vegetation, soil 
moisture, or soil particle size and soil structure. Examination of the 
Zeiss and AMPS imagery (see figures 22-39) reveals no vegetation in 
field D. Differences in soil moisture between field D and the other 
fields is possible but somewhat unlikely since (1) no measurable rain 
was recorded during the three days prior to the aircraft overflight, and 
(2) there is no field irrigation in this area during the late fall 
(e. g., October) . During this time, however, many of the fields undergo 
tillage treatment. Stoner and Horvath (1971) have demonstrated that 
soil surfaces disturbed by tillage treatments exhibit a decrease in 
reflectance for the rougher soil surfaces. Other studies (Bowers and 
Hanks, 1965; Orlov, 1966; and Shockley et al., 1962) conclude that j 
increasing particle diameter results in a decrease of reflectivity. 
Therefore the statistically higher radiance values in field D could be 


B-93 


n**s CM*k«ii « tlB( t 1 lAM % 

•I •» It W 4 i i» *1 M •» •« *• «I^A •• '' *• •• •* 

i« I* i» I* i»\*i •• *' *• •* •' *K?* ** II n i» H ^ 

•» II I* !• /♦ l» l» ** ** ** *' ** ** II Q* y !• I* !• •! II 

It I* i» r*!i* t*^i^ i» i\^ii mS\ fiy •» »i *1 •! •• *♦ •» II !• M i» y •• 

It It It tjt* /t(^<* #»\# || •< •» ♦» •» •' •• M It II y •» ** ♦v'l It 


II It II II i» n n ri »t i\y I 

It n^l iKi* !• It l« It It 1 

II vJ* iiV II r* ft /• f* 

i« n II It mit ft It |t |» lO 


lK«t I* tl «l •! tt tl 4 


•t tif A| tl tl It II 
*t ti^ t* I* 


II II II i\ I* It It y •• 
^ It II */ If It f*A% It 


i^t It t\ii II II II II 
I ^ It II Xn II II II 




II il t« II II It II II II |(||p 

II o>it It II il II II II l^.in 

II II jV It II iiV^ II II 


II It k% II II II II II It 1 ^ 

'T~y^ i*\ *1 II It It II i| 



»l II II II 1 ^ It II II II II It 1 / 4 
II II 11 It Il II II II II II j 
II II tf II W OT II II It II It II 4 
II II II II A y II !• Il II II It I 
II II II II II It II II It II • 


l» Tt 1-^ II %\ * 

III I. (Trrx^ 




’» ’» •» •» 

•• iimii II II II II 


Mr ** ^ ** ******** *f 

:'^ki It It II It It H •• •» •» 4 

I ^iki II It II II II It II 1*0 II II II iK, 


>1 11 II It II It It II II |l II U U 4 


lt>% II tl I II II II k lllkl II 11 II It II II 11 /^ II II It II II I I I I I 

It iijii Y^'*\ ** ** n^li * ** ** ** *' *^ **l** V ** ** *• jl* ~ \ * * * 

It iilii I I V I II •' *♦ ** •» •* ** *»W •» 

It Liljt ,i,t >w, ti.L'i.v It It It II II II II X'*-* M It lya^mX* » It I 4 

l••lV^ v*/i»i 4 * V * ** *' ** *' /*• 

li It )» ll\ » N* • w •’ •• ** •' •» •* •* H **V* ** •» •* Wh ••/*'* 

'V'V * r* •’ •' *' •’ *’ ** *'/'* *• V*’ ** ** l|i^* It rt y I II 

II \ tt\»l \ I I I t^lkl^lt i< II II 11 i^l* II I^VM *' II •• •tio/^fl 

II nVu^iy ii\ t I I i[4!Lt i«ii II li 11 II II y It It 1^ Il ii{yyrr-»v\i f I 


II II t\ ityi \ I I I tM,ki^ii i< II II II i»\i* II iiNm *' II 

It II ir^t^iiX til if4!»i itif II II li 11 11 i>y It It if^ Il 

nts^ 

II II •Y^**|' * * * ipj^'r*!*' •* *' *’ *’r* '* ** •’ 

Wt/ II tl nvi*'ii y I I I j||)Ji Il u •• *1 I'j 1 * i» !♦ ii 


1 i 6 !v '* *' tl II ii I */^ jf'r p* y *' *y^ *’ *1 *• '• *' *' '• Y 

*IVii II II II II ii\i *1 r If •* I' »• I' **yt 


ty^\ II w II II 
''^1 II U 11 tl II 


II II «i » 
II II iJ t 


> ' ’ 

(f I iVi> I t 



I t f I f I I 


W Y\ I I I 
»i A\t I 1 I 


ii^t It tt\iA f II 

9 \ t» If tf • • tl 11 


II tl/t It t* I I 

I t* II *M’'\ » » 


n il •! »l 1 ^ «« » l I* It 

ii i« •* J It !• 

ly ly At It %f *1 X* !• If 

kl /t M It >t,Vi k 1 It It It 


• It If y t'^t It If It It |4 »i $• *1 

> It >t^ •/it It It t« It /i »t ft *t Ui^l 

I If If Nt ti\it It it i» !• iSk* »y\i tl It II Ui Ux» 

• It It i*\i X »• '• »• *' '• y t^ )*’ If K tl 

» It It !• 4 tij i« If It !• It 4 •• #iiyi/*i $f »y •• 

r\> It It iA*y It If It It iy»i t If 
V* It It it^ It It It It iA*i • Pyti y 

' It It It It ^ It If I# It »t\i «• mt 1^11 

» It It I* It It It It It If If y I 

»• *' *• It It If If It It yi A ** j* * 

h •: ’! ” ” !! !: :: :: ::/• rw .y^ <*iH 


} It It It It It It It It It It 4 . A •• < 



nh^ » 4 1 *| It It It It It It II It If I* ) 

|ljA 11^ It If If It It It If If It !• 4 

^(;|U If It It If If •• It If If If |» ttU 

If It If If If It It II II It M 

j^jilYt M It It It It It If If If |f 

r1^ It If w It It If II It If i4 « 


Jiy It If w It It If II It If I* k« f ^ 

'riY Mo** *' *' '• *• •* *' *' '^<r fl'W ** 

. •• •^I'^t It It If It If I* If If It yi •• tf tf 

uiy It It It It i« It It If It if i^t* it^i • tt tl 

p,'i\ tl t# It •• It It It If If If If W »i ^viWf 

pl^Mki tvt It It It It if It If It iJ •• » f 


( tl t« t\lt It It !• |f II If If 


It If t« It If If If |4 I 


n •••i 

!fl »f M tl tl •« tl t<\l> It It ^ 



w! »i tl tf *1 tl tl tl t«/it It It if^ 4 
W 11 tl #1 tl tl tl tl •'( It If If It A 4 

k|jU tl tl *^ i» It u If i»j< 

e;» tl tl If If If If If y ( 


I || tl At It It It It Ai 
r^ly^lt It It If If \t 
I It •• If If If It If l\ 


l»i tl tl ^It It •« It If If II ltA» 

ki tl tl ti^ Il It It If It It If V 

iy tl tl *• X** »' '* *' ** *' *' y 

lAkf •! tl tl tl^M^ It If If If /♦ 


I 1 I I 

Inz 


Contoured radiance values for the channel 9 (0. 760-0. 860iira) 
of the nultispectral scanner. Field B is shown. 










TABLE 17. M^S Spectral Band vs. Field Radiance 


Average Field Radiance (Watts/cm^-Sr-AX) 


Spectral Band 
(urn) 

Field A 

Field B 

Field C 

Field D 

X 

S 

.440-. 490 

130 

133 

145 

145 

137.8 

7.4 

.495-. 535 

120 

121 

130 

135 

126.5 

7.2 

.540-. 580 

93 

95 

102 

107 

99.3 

6.5 

.620-. 660 

83 

83 

90 

96 

88.0 

6.3 

.660-. 700 

82 

82 

92 

100 

89.0 

8.7 

.760-. 860 

200 

205 

225 

252 

220.5 

23.6 

.970-1.060 

137 

152 

160 

175 

156.0 

15.9 


accounted for by recent tillage in the other three fields. However, as 
mentioned previously, no ground data is available to verify this 
postulation. 

A comparison of the M^S-derived radiance vs. wavelength plots of 
fields A, B, C, and D (see figures 42-45) to Condit's laboratory curve 
of reflectance vs. wavelength for chernozem type soil (see figure 40) 
reveals several differences. The MS data shows a decrease in radiance 
from 0. 440}im-0. 700ym, an increase from around 0.700)am to O.S60pm, and then 
another decrease from approximately 0.860pra to l.ObOnm. This contrasts 
v/ith Condit's curve which shows a general increase in reflectance vs. 
wavelength from O.SOOjjm to 1.000|jm. The differences are not totally 
unexpected as field spectra often have differences with reflectance 
spectra from soils measured with spectrophotometei’ measurements in the 
laboratory. These differences arise because of (a) 02> CO 2 , and H 2 O 
absorption reducing incoming solar radiation in certain spectral bands, 

(b) variable solar illumination in different conditions, (c) soil 
texture and structure, and (d) solar intensity peaks at O.Sym, falling 
off rapidly at shorter and longer wavelengths. Even though the labora- 
tory and field curves do not correlate, it is encouraging to note that 
the M^S data provided an identical spectral "signature" for four fields 
of the same soil type. 

C . M^S Thermal IR Data and PRT-5 Data 

Several experiments (Blanchard, et al., 1974; Idso et al., 1975; 
and Reginato et al . , 1976) have been performed demonstrating good rela- 
tionships between soil-surface temperature and soil-moisture content. 

A summary of these reiationships is found in Schraugge (1976b) . Since no 
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ground truth was available for the October 17 experiment, no comparison 
was possible between it, the PRT-5 radiation thermometer, and the M^S 
channel 11 thermal IR data. Because the "footprint" of the PRT-5 sensor 
is approximately .11 km and line 8 of the aircraft data was not centered 
over the test line, the PRT-5 data does not include any of the fields 
A, B, C, or D. Therefore no comparison is possible bet»\?een the PRT-5 
temperature data and the temperatures derived from the M^S thermal IR data. 
D. Landsat Data 

The closest date for useable Landsat data for the October 17th 
experiment was October 1st. Since the test date and Landsat overpass date 
were not coincidental, the only comparison which can be made between the 
two is to determine vvhether the individual fields A, B, C, and D, which 
are visible on the M^S imagery can be distinguished on Landsat imagery. 

The purpose of the comparison is to determine if Landsat data can be used 
as an aid in studying individual fields in the Luverne, Minnesota test site 

Figure 46 is a side-by-side comparison of the channel 4 band (.54- 
2 

.58iim) of the M S data band 5 (.S-.6pra) of Landsat. The Landsat data 

* ' • ■ ■ ■■ 9 ■' ' ' 

has been photographically enlarged to approximate the scale of the MS 
aircraft imagery. Field A through D have been labeled on the M^S image 
as well as other scene location features . The Landsat imagery has 
been similarly labeled. The result is poor field differentiation in 
band 5 for this area during October. None of the test fields can be 
located. — 

Figure 47 is a side-by-side comparison of the channel 10 band (0.97 
~1.06vim) of the M^S data and band 7 (0.80-l.lym) of Landsat. Once again 
the test fields and SGene location features have been labeled. Even 
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LANDSAT(.5-.6mn) 


A comparison of the channel 4 band (0.54-0.58um) 
band 5 (0.5-0.6ym) of Landsat-2. Identical feati 
bands are numbered. I 









Fig. 47. 


LANDSAT (.8-l.lum) 


A comparison of the channel 10 band (0.97-1. 06pm) _of_the.M^S 
and band 7 (0.8-1. 1pm) of Lan dsat-2. Identical features on 
both bands are numbered. ! 


througii the test fields are not shov^’n in band 7 of Landsat, it is much 
easier to distinguish individual fields during October in band 7 than in 
band 5. The question still arises as to why the test fields A through 
D do not appear on the October 1st satellite imagery. Other bare fields 
(e.g., 5 aJid 9) the size of the test fields do appear in band 7. According 
to the climatological data for October (NOAA, 1975) there was only one day 
of appreciable rainfall (.21 inches on October 9th) batween the time of 
the satellite imagery and the subsequent aircraft imagery. Thus, different 
moisture conditions cannot account for the difference. The answer may be 
because the fields had not yet been tilled by October 1st after the fall 
harvest but were by October ISth, the time of the aircraft data. 

Conclusions 

Owing to the lack of ground data during the October 17, 197S experi- 
ment, no conclusions can be drawn concerning soil moisture measurements 
using remote sensing techniques. Three conclusions, however, were 
derived from the remote sensing data. They az'e (1) multispectral data 
can be used to assi,gn a spectral "signature" to specific soil types, 

(2) knowing prior meteorological conditions, statistical analysis of 
multispectral data can be used in nonirrigated areas to give a qualitative 
determination of field roughness, and (3) Landsat MSS band 7 data can be 
used as an aid in remote sensing studies of individual fields which are 
within the resolution limitations of the satellite data. 
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